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Abstract 

An example of application of the specialized computer algebra system GRGec to 
the searching for solutions to the source-free Maxwell and Einstein-Maxwell equa- 
tions is demonstrated. The solution involving five arbitrary functions of two variables 



(2jT), is presented in explicit form (up to quadratures) 

> : 

X : 1 Introduction 

The specialized computer algebra system GRGec [1] is intended for applications to the 
classical theory of gravity as well as adjacent problems related to the classical theory of 
field and geometry. In particular, GRGec is currently 'aware' of the majority of main 
characteristics of the geometry of a curved space utilized in Einstein' theory. These are, 
for example, the bases in foliations of exterior forms connected with metric, the connection 
(including its Newman-Penrose representation), the curvature, its irreducible constituents 
and algebraic invariants, the general equations connecting these objects (Cartan' structural 
equations, Bianchi' equations, various algebraic identities), the field equations of the grav- 
ity theory (Einstein equations for vacuum and various matter content). This list can be 
continued with the basic elements of the Rainich theory, the theory of Lanczos potential, 
the methods of the symmetry description (Killing equations). GRGec system operates with 
major characteristics of such classical fields as electromagnetic field, massless spinor field 
(Weyl field), massive spinor fields (Dirac fields) including the both cases of the inclusion 
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of interaction with electromagnetic field ("charged" Dirac field) and the absence of elec- 
tromagnetic interaction, massless scalar field, conformally invariant scalar field, massive 
scalar field, massive vector field (Proca field), pressure-free dust matter, massive and null, 
etc. 1 

Concerning the above characteristic of GRGec capabilities, it has to be noted that there 
is, of course, a number of programming packages of symbolic manipulations with a similar 
application field (an excellent survey of the latter topic is given in Ref. [2], see also [3]). At 
the same time a general overview of the current state of this field of computer applications 
led a majority of authors to the conclusion that, in spite of noteworthy achievements of 
the applied computer algebra, it would be still premature now to single out any of the 
existing programming systems of the class implied. None of them can be assessed as an 
universal tool which reveals unconditional advantages over all the potential competitors and 
best suits for majority of applications. Accordingly, in practice, answering what is the best 
system of symbolic calculations for application to the gravity theory, one has to preliminary 
make definite what a class of problems to be tackled is implied. This observation allows 
one to suppose that GRGec can find an own position in the row of computer-based tools 
utilized for theoretical investigations in the field of the gravity theory. 

It is worth mentioning that one may find a lot of examples giving evidence that the 
practical usefulness of programming instruments intended to 'non-local' applications is a 
fairly subtle matter affected by a plenty of factors, often rather vague and sometimes fairly 
unexpected. In particular, it may be stated that any partial characteristic (similar to 
the above list the physical-mathematical material implemented in GRGec) or any verbal 
description is not able to ultimately establish a proper measure of the actual advantages 
— as well as drawbacks — of a software given. 

Nevertheless, there is an evident straightforward mean to gain some insight into the 
heart of the problem. Specifically, an overview of a proper set of manifest examples of the 
practical use of a programming system can serve a probe exhibiting its real characteristics. 
Concerning GRGec system, the present work realizes a step just in that direction. 

In order to demonstrate in acting the basic features of GRGec system (and to exhibit 
some of its capabilities) we would like to learn the solving of a concrete problem, namely, 
the integrating of the system of Maxwell and Einstein-Maxwell equations. The result 
we shall derive is of a certain independent interest. The solution of the electrovac field 
equations obtained below with the help of — or, if one prefer, in collaboration with - 
GRGec system involves five arbitrary functions of two variables. Though the calculation 
we picked does not pretend to be related, with regard to its complexity, to the category 
of record ones, it nevertheless gives, in our opinion, a clear evidence of a usefulness of the 
programming tool presented for the handling of such sort problems. 

It is worth noting that our choice of the subject of consideration was severely condi- 
tioned by fairly simple 'pragmatic' reasons. On the one hand it is evident that consideration 
of a toy problem would obviously decrease a potential interest of the demonstration. In 
any case, it tells a reader nothing definite in connection with the real efficiency of the 

1 This list should not regarded as ultimate one. 
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computer tool applied. On the other hand, in the opposite extreme case of complicated 
problems lying on the edge of the practical capabilities of the method, there would appear 
an evident technical obstacle strongly restricting the scope of discussion. Indeed, solutions 
of complicated problems are usually too bulky to be exhibited in sufficient details within 
frames of an article of any reasonable length 2 . In the best case the only initial posing of the 
problem and the final result with minimal discussion can be presented, the explanations 
of the essence and, all the more, the features of the course of the solving procedure being 
inevitably omitted. (As an indirect evidence of such a state of things we might refer to a 
lot of publications of investigations resorting to assistance of the computer algebra where 
a remark "... with the help of computer algebra we found that. . . " , or the like, is the only 
'discussion' of the relevant programming issues.) 

In our case, as usually, a contradiction of the choice among two opposite extremes is to 
be settled by means of a compromise. Specifically, although the problem we shall discuss 
below is in principle tractable by means of a 'manual' calculation (being rather error prone 
though), definitely, nobody would qualify it at whole as a trivial one. It seems thus to 
be able to provide a proper material for the displaying the fashion of the work with the 
computer tool considered. 

Reverting to the general characterizing of GRGec system, it has to be mentioned that 
one of the major goals pursued during its development was an attempt to release a user 
from the duty to 'develop a program' — at least, in the meaning usually associated with the 
latter pursuit. Rather, dealing with GRGec system, one has to merely describe, employing 
mostly 'habitual' words arranged in a sequence of natural 'phrases', and operating with 
more or less standard mathematical notations, what initial data are given and indicate what 
a result has to be generated. To that end, GRGec maintains the input language which is, 
probably, as close as possible to the one used for the representation of the relevant notions 
and the relationships taking place in the application field itself. Of course, comparatively 
sophisticated mathematical manipulations require from a user an additional control over 
the calculations. In particular, only he or she is able to decide what a way should be most 
fruitful under the specific conditions depending on the details of the method of the problem 
treatment and the data given. However this, anyway, means that a researcher may focus on 
the essence of the problem considered, the technical programming-related issues requiring 
considerably less outlay than the other symbolic computational tools would necessitate. 

Currently, GRGec system enables one to calculate or to subject to the other processing 
more than two hundreds of the so called data objects modeling the basic notions and 
relationships (equations) originated from the field theory in the curved space-time and 
the geometry. However GRGec is n °t accommodated for the purpose of the abstract 
index manipulating. The handling of data objects endowed with extrinsic indices (tensors, 
spinors, etc.) is carried out utilizing, essentially, the explicit sets of their components 
attributed to a definite gauge. 

Touching in brief the topic of realization, it is to be noted that GRGec is a 'superstruc- 

2 It may not be definitely stated that, after all, this characteristic of the present paper proves to be 
unconditionally satisfactory but there is seen no ways to amend it in that respect. 
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ture' over the well known general purposes computer algebra system Reduce developed 
by A. Hearn (see Ref. [4]). Accordingly, GRGec potentially possesses the same degree of 
portability as Reduce does. The programming language utilized for the GRGec realization 
is not Rlisp, the language of the own Reduce coding, but the Lisp dialect named standard 
lisp, see Ref. [5], which is supported in frames of the PSL (Portable Standard Lisp) pack- 
age, the bottom level of the Reduce infrastructure hierarchy 3 . It is important to emphasize 
that a user need be familiar with neither standard lisp, nor with Reduce (although a 
knowledge of the latter system would not be completely useless of course especially in the 
cases of more sophisticated calculations). Specifically, GRGec creates an 'opaque casing' 
which mostly conceals the lower level data structures utilized and the inner ways of their 
handling. Reduce is thus a 'lower level' structure from viewpoint of GRGec- 

A language exploited for the communication of a user with GRGec system is a product 
of the independent working up. It is closely related to the application field and, at the same 
time, mimics the elements of the natural language, utilizing a sort of partial simulation 
of English. Accordingly, the essence of GRGec programs is, as a rule, understandable for 
any specialist in the application field — even a dilettante in the programming as such. In 
the worst case, a quite moderate number of elucidations of specific issues could prove to 
be necessary (one will see this on examples considered in the present work). 

The source code of GRGec amounts now to approximately 1.5 Mb. The compiled (32- 
bit) binary code occupies about 1 Mb of a disk space. GRGec is currently intended for a 
free distributing. 

2 Preliminaries 

It has to be noted that in the case of four dimensional Lorentz geometry GRGec uses 
for calculations a version of the null tetrad method based on Cartan' structural equations 
incorporated with 2-spinors theory. 

Unfortunately, it is problematic to point out a single source which would survey the 
formalism coinciding in all substantial details with the one implemented in GRGec system 4 . 
At the same time the univocal representation of a mathematical matter underlaid is of 
a notable importance for anybody who deals with a concrete computer realization of a 
mathematical formalism. 

Specifically, a well known detailed presentation (now, essentially, unofficial standard) 
of the methods of spinor calculus is given in the well known monograph by Rindler & 
Penrose [6]. However, at first, the methods based on Cartan' structural equations and 
their incorporating with spinor approach are not considered there in fact. At second, apart 
from some notational discrepancies, there is a distinction in the conventions fixing the sign 
of the spinor contraction, i.e. the definition what a rule, either la = £ab t B or la = l b £ba, 
is to be used for spinor indices manipulations. The choices adopted in [6] and in GRGec 
are opposite. [The point is that the both above definitions were (and are, see, e.g., [7, 

3 It is worth noting that PSL is not the exclusive possible base suitable for the Reduce implementation. 
4 This does not mean of course that the latter is unique in any respect. 
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Appendix A: A]) utilized in the literature. The former one is adopted, in particular, in 
the series of investigations by J. Plebahski et al considering complex Einstein spaces, c/., 
e.g., [8], Eq. (3.22).], this system of notations and basic conventions being elaborated a 
decade earlier than Ref. [6], adopting alternative choices, appears 5 .) Further, the methods 
of the calculus of exterior forms in a version adapted to the gravity theory are outlined in 
Ref. [9]. However, the spinor based description of the geometry and field theory is given 
there after a fashion fairly incomplete and insufficient for our purposes. We might also 
refer to the work [10] but it treats a specific problem and cannot serve a review of the 
formalism applied; besides, there are still some disagreements in conventions utilized in 
Ref. [10] and GRGec, unfortunately. 

We have to give therefore here a summary of the main elements of the mathematical 
methods which will be 'applied by GRGec system' for the treatment of the problem con- 
sidered in the main body of the paper. Besides, we discuss here its mathematical posing 
necessary for the subsequent processing by the methods of computer algebra. 



2.1 Elements of mathematical formalism 

Let 4-dimensional Lorentz metric be represented in the form of expansion 

Here and below undotted spinor indices (denoted by upper-case Latin A,B, etc.) run over 
the two element set {0, 1} while the dotted ones (dotted upper case A, B, etc.) run over the 
(distinct) set {6, 1} of 'dotted integers' (subjected, of course, the standard arithmetic). We 
use throughout the standard spinor notations normalized in accordance with the following 
convention which determines the spinor index manipulation rules [8]: 

/ — c i B i B — i C AB i — c , B , B — , A B 

b A — t AB b i h — i h A — AB ' — A 

Here 2-dimensional Levi-Civita symbols e**, e** may be defined as follows: 

c — c c — 1 • C AB — C BA c 01 — 1 ■ 

£ab — —£ba, e i — 1, e — — e , e — 1, 



~AB ~ BAi c 01 



^B = _ e BA e 01 = L (2) 



Further, i) A ^ (involved in Eq. (1) with distinct positions of indices) is the tetrad of the, 
generally speaking, complex- valued covectors (1-forms) which are null with regard to the 
metric g. In the case of the metric signature ( — h ++) assumed throughout they always 
can be picked on to satisfy the Hermitean-like symmetry relation 



»0d=*00> »H=»1V »0i=»10- ( 3 ) 



5 In connection with these remarks it should be mentioned that the predecessor of GRG EC is dated to 
1981-2. 
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In such a case, the two tetrad elements ^QQi^^j are rea l while i?^ , i? ^ are necessarily 
complex. (The meaning of distinction in the format of numerical indexing which is observed 
in Eqs. (2) and (3) will be explained below.) 

We shall refer to Eq. (3) as the real gauge condition. The imposing of this gauge does 
not completely fixes the tetrad: it may be still undergone gauge transformations whose set 
includes, in particular, the proper orthochronous Lorentz group. 

Further, it is sometimes convenient to utilize another indexing of the tetrad. Let us 
consider the following re-definition 6 : 

& = *ov §l = § id> & = *od> ^ 3 = -^i- ( 4 ) 

In terms of such an indexing Eq. (1) takes the form 

g = 2tf° (s)^ 1 + 2tf 2 ©tf 3 , (5) 

where a(§)/3 = \ol <g> j3 + \f3 ® a denotes the symmetrizied tensor product. We shall say 
that $ A] j is the tetrad endowed with spinor indices while Eq. (4) defines the tensorial (or 
tensor) tetrad indexing. 

Important algebraic objects which characterize the property of space-time metricity 
are the two families of the so called S'-forms constituting 2-index spinors. In particular, 
undotted S-forms are defined as follows: 

Sab = \ e k H Ak N® BL = \ 9/ A fl . . (6) 

Next, each metric induces a unique torsion- free metric compatible connection. Within 
framework assumed, it is described by the object endowed with a pair of symmetric undot- 
ted indices (which do not constitute a true spinor however) whose components are 1-forms. 
It is named undotted connection and is denoted Tab- Due to indicial symmetry there are 
maximally three distinct components of undotted connection. They are denoted T q, T T «?. 

Notice that we used above so called summed spinor indexing which applies to any 
object symmetric with respect to all spinor indices of the same class (undotted and/or 
dotted ones) including the case of connection. Hereinafter, dealing with numeric values of 
spinor indices in explicit (expanded) formulae, the summed indices are used throughout 
(see Eqs. (3) and all the subsequent ones; Eqs. (2) are an exception). In particular, in 
the case of the undotted connection the complete collection of its components comprises 
r o = Tq,Tqi = Tio = Tg, Tn = T2 (here at right hand side the summed indexing is 
used while at left the corresponding 'numerical values' of 'ordinary' spinor indices are 
substituted, cf. Eqs. (2)) 7 . 

Given null tetrad, the connection forms can be computed from the so called first struc- 
tural equations which are also named the first Cartan equations. There exist two their 

6 Within framework of the matrix language this point assumes introduction of Pauli matrices. We do 
not utilize them in an explicit form. 

7 As far as we know the simple but very useful notion of summed spinor indices was never explicitly 
introduced and used on a regular ground. 
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representations. The first version of the first structural equations involves the differentials 
of the tetrad elements and, in the case of the fulfillment of the real gauge condition, reads 

d^ + r$A^ B +ff A9 AC = 0. (7) 

Here the complex conjugation varies the class of connection indices, i.e. enters the 
equation, effectively, with two dotted indices, B and C, which are connected with undotted 
B and C, respectively, by means of the map sending to 6 , 1 to i and vice versa. [More 
precisely, the so called dotted connection forms T^, which are in general case algebraically 
independent on undotted ones, should be situated in Eq. (7) instead of T^. However, under 
the real gauge condition, and T B are mutually 'Hermitean conjugated' and we prefer to 
evade introduction of an additional object which is, essentially, superfluous for our current 
purposes.] 

The second version of the first structural equations (first Cartan equations) involves 
differentials of S-forms instead of ones of tetrad elements and reads: 

&S AB + 2T<( A /\S B)C = Q. (8) 

The both versions of the first structural equations (which are, essentially, equivalent) always 
can be algebraically resolved with respect to the connection forms. 

The next important geometric relationship utilizes the so called undotted curvature 
spinor VL AB — ^ba whose each component is a 2- form. The link of undotted curvature 
and undotted connection is yielded by the second structural equation 8 . It reads 

n B A =dT B A + TfAT B K . (9) 

Space-time is locally flat if and only if Vt A (and its dotted counterpart) vanishes. 

The specific relationships lying in foundation of the general relativity (Einstein equa- 
tions for the general matter content and their particular case, Einstein-Maxwell equations, 
which we shall deal with in the present work alone) are formulated making use of the so 
called irreducible spinor curvatures. The latter involve in particular, the undotted Weyl 
spinor ^abcd = ^(abcd) representing Weyl curvature tensor and its Hodge dual, Ricci 
spinor <& ABCD = &(ab)(cd) equivalent to the traceless part of the Ricci curvature tensor, 
and the scalar curvature R. All these objects can be determined from the known undotted 
curvature (and S-forms) with the help of the equations displayed below: 

i ^abcd ■ Vol = Vt (AB A S CD ), -i ®abcd ' Vo1 = ^ab a s cd> ( 10 ) 
iR-Vol = AVL AB AS AB . 

Here the dotted S'-forms S^^ are defined similarly to undotted ones by the equation 

Sab = \ eKL ^KA A $ lb 

8 undotted one, there is also dotted second structural equation involving dotted curvature and dotted 
connection 



7 



(cf. Eq. (6)) and 



Vol = i- 2 9 B A A6 c B A6°A6A 



is the volume element (nonzero 4-form) induced by the corresponding metric. 

Now let us outline the form of the basic equations of electromagnetic theory within 
framework of the formalism under consideration. They are extremely transparent. Elec- 
tromagnetic field is associated with the complex- valued 2-form uj which is named undotted 
electromagnetic 2-form. It has to be spanned by the undotted S-forms, i.e., admits the 
expansion 

cu = cj> AB S AB . (11) 

The above coefficients <j> AB = <p BA constitute the so called undotted spinor of electro- 
magnetic field. The charge- and current-free Maxwell equations reduce to the closeness 
condition imposed on the undotted electromagnetic 2-form: 

du = 0. (12) 

We shall refer to it as undotted Maxwell equations. In the case of the imposing of the real 
gauge condition (3) the complex conjugation of the undotted electromagnetic 2-form yields 
the dotted electromagnetic 2-form uj = uj which admits a similar expansion 

Co = cj) AB S AB , 

where the coefficients <p AB constitute the dotted spinor of electromagnetic field. Finally, 
under a suitable choice of the physical units, Einstein-Maxwell equations are represented 
as follows: 



2 ®abcd = 4> AB 4>cbi ( 13 ) 
R = 0. (14) 

Eq. (13) is called the spinor part of Einstein-Maxwell equations while (14) is their scalar 
part. 

Now we may proceed with the description of the problem chosen for the demonstration 
of the application of the formalism outlined with the help of GRGec- 



2.2 Radiative electrovacs 

In the present work, we shall consider solutions of the source-free Maxwell and Einstein- 
Maxwell equations (electrovacs) which, under a proper gauge, satisfy the equation 

T = 0. (15) 

Characterizing its geometric, the first important consequence follows from '-component 
of Eq. (7) which now reads 

d # 2 + 2 # 2 A Re = 0. (16) 



8 



Specifically, by virtue of Frobenius' theorem, 9 2 determines an integrable distribution. Said 
another way, there exists, at least locally, a scalar function x such that 

6 2 = E l( ix (17) 

for some non-zero factor Si, the hypersurfaces x =constant being null. 

Further, it has to be noted that the coefficients p, a, . . . ix of the general expansion of 
undotted connection with respect to null tetrad 

T = p0° + oQ 1 + rQ 2 - k6 3 , 

T i = a6° + (36 1 - 7 2 + s9 3 , (18) 
T 2 = \6° + (10 1 - vQ 2 + tt6 3 . 

are nothing else but the Newman- Penrose (NP) spin coefficients, see Ref. [11], [9]. Some 
of them (namely those involved in expansion of Fq and Fj) describe the invariant prop- 
erties of the null congruence determined by the covector 6 2 , i.e., generating hypersurfaces 
x=constant. In particular, if k=0 (that follows from (15)) this congruence is geodetic and 
possesses the complex expansion equal to p and shear \a\. Thus we immediately see that 
all the optical scalars vanish in the case under consideration. 

Of course the ansatz (15) imposes severe restrictions on the space-time curvature as 
well. Specifically, the '^'-component of Eq. (9) reads 

n = dr + 2r Ar i (19) 

and, thus, Eq. (15) immediately yields Qq = 0. This, in turn, implies in particular the 
following constraints on the components of the Weyl spinor and the scalar curvature (see 
Eq. (10)): 

= m = ^ 1 = m 2 + ±R. (20) 

Since in electrovac space-times the scalar curvature R vanishes (Eq. (14)), \l> £ vanishes as 
well and only ^ j may be non-zero. 

The above restriction on the components of Weyl spinor implies that the space-times in 
question belong to Petrov-Penrose types III or N (or, perhaps, are conformally flat, the type 
0). The corresponding multiple principle null direction coincides with the one determined 
by the congruence of null geodesies generated by null hypersurfaces x =constant. We have 
seen that all its optical scalars vanish. 

It is known that gravitational fields with the above algebraic structure of the conformal 
curvature are usually associated with gravitational radiation spread in the 4-dimensional 
light-like direction, multiple principal direction. In particular, the parallelly propagated 
plane waves (pp- waves, see Ref. [9], section 21.5) satisfy all the above conditions. Due 
to these reasons and for convenience of further references, we shall name the solutions 
of the source-free Maxwell and Einstein-Maxwell equations satisfying the restriction (15) 
radiative electrovacs 9 . 

9 It should be mentioned however that there exist configurations describing radiation which do not 
belong to the class singled out. A wcllknown example is, in particular, Robinson- Trautman' class of 
metrics, see Ref. [9]. On the other hand, there are 'radiative' spacetimes (in our terminology) which are 
unlikely to be associated with any radiation process, e.g., the Hauser' solution [12]. 
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In the present work, we find the 'most general' metric belonging to the class of radiative 
electrovacs, reducing all the field equations to several subsequent quadratures. Our solution 
does not actually includes all the radiative electrovacs but the metrics missed here (e.g. 
pp-waves) are either generated from it by means of appropriate limiting processes or arise 
following a particular 'degenerated' side-branches of the 'generic' integration procedure. 

Now we are able to carry out a complete posing of the problem. 

2.3 Posing of the problem 

First of all we have to determine what a tetrad fits the restriction (15). To that end, we shall 
utilize some further implications inferred by Eq. (7). Specifically, it has been mentioned 
that its '^'-component implies Eq. (17). Absorbing the factor Si by the tetrad element 
9 3 (which is an admissible gauge transformation, cf. Eq. (5)), we obtain that, without loss 
of generality, one may assume 

9 2 = dx/V2. (21) 
Further, eliciting the '^'-component of Eq. (8), one obtains 

dSg + 2S AT 1 = 0, 

where Sg = 9 2 A 9° (see Eqs. (4), (6)). Again, the Frobenius' theorem implies 

Sg = e a dx Ad0 

for some functions a, x, 0. By virtue of (21), without loss of generality, one may get x = 2 x 
(we a priori know that by virtue of definition Sg = 9 2 A . . . = (...)dxA...) and, then, 
executing, if necessary, some gauge transformation, we obtain 

9° = e a d<f)/V2 9 1 = ¥ = e«dj)/y/2. (22) 

It is immediately clear from Eq. (5) that the only value of (Re a) is significant while the 
choice of (Ima) is a matter of gauge. Indeed (Ima) may be arbitrary, leaving the metric 
unaffected. Additionally, analyzing the above explicit equations, it is easy to see that the 
local coordinate transformation — > E 2 ((f),x) for arbitrary S 2 holomorphic with respect to 
(and satisfying the natural restriction S 2 ,</,7^ 10 ) still remains allowable, provided the 
corresponding transformation of a accompanies it. We reserve this possibility for a future 
use. 

Noting that 9° A9 1 A9 2 ^ 0, the variables x, 0, are functionally independent (and thus 
is complex and 'functionally independent' of 0, i.e. and are a form of representation 
of two real independent coordinates). Let us use them as three of four coordinates and 
denote the fourth independent real coordinate y. 

Thus, up to now, the only non-specified tetrad element is 9 3 . Its generic holonomic 
expansion reads 

y f 29 3 = S 3 d?/ + 2Re(xd0) + hdx, 
10 Hereinafter, a subscript following comma denotes partial derivative. 
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where the coefficient H 3 may not vanish, x (complex) and h (real) are some functions. 
Then, executing the admissible change of the coordinate y, y — > f S 3 dy, one always may 
reduce (the counterpart of) S 3 in new coordinates to a constant. Thus we may assume 

^20 3 = dy + 2Re(xd0) + hdx. (23) 
Let us notice further that the transformation 

y -> y + 2 4 (24) 
for arbitrary real function H 4 = H 4 (x, 0, 0) which is accompanied with the replacing 

X^X-2 4 ,^, (25) 

and some transformation of h is also a change of a gauge. (It does not affect any of 
equations considered above including the expansion (23).) Notice also that then we have 
the following subsidiary transformation; 

X,4,+X,<f^ (x,$+Z<i>) ~ 2S 4,^- (26) 

It will be useful below. 

Ultimately, we can formulate the problem as follows: We shall search for electrovacs 
whose metric is described by the null tetrad (22), (23) and possesses connection which 
satisfies the constraint (15). 



3 Handling field equations with the help of computer 
algebra 

In this section subdivided into several subsections we present the application of GRGec 
system for the reduction of the field equations specifying radiative electrovacs. In each step, 
a separate series of calculations is carried out yielding some new information describing 
the structure of solution. Whenever its 'amount' exceeds some reasonable threshold, the 
calculation stops and further proceeds with the next step from the beginning, utilizing the 
relationships obtained so far as a new and more detailed initial data. 

Any knowledge of GRGec system is not presupposed. Accordingly, relevant expla- 
nations on this issue are given whenever necessary 11 , making the discussion, essentially, 
self-contained. 

3.1 Step 1: Coding of the initial data and inference of first im- 
plications 

Generally speaking, the data to be read and processed by GRGec is to be arranged in 
a form of the so called problem specification which is a regular text usually stored in a 

11 excluding some general means of the formulae coding which are utilized in majority of programming 
languages such as Fortran, Pascal, etc. 
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disk file. (One finds in Appendix A a sample of the workable problem specification.) The 
structure of problem specification is fairly transparent. Apart from a title and conclusion, 
it consists of a series of sections comprised, in turn, of several paragraphs and must, in 
particular, contain the section of (initial) DATA. The other sections (for example, the section 
of INSTRUCTIONS, the section of SUBSTITUTIONS, and some others) are optional. 

Specifying initial DATA, a user has to explicitly list all the non-implemented ('user- 
defined') objects which he or she is going to utilize. Here the particular symbols of C00R4- 
DINATES, specific FUNCTIONS ('unknown' to the system), etc., are, in particular, meant 12 . 
In our case all what we need at the beginning is explicitly DECLAREd and defined as follows: 

1 input: Data : 

input: declare COORDINATES x,y ,phi ,phi~ ; 

3 input: declare SCALARS al , chi ,h(x ,y ,phi ,phi~) ; 

input: declare REAL x,y,h; 

5 input: declare COMPLEX CONJUGATED phi & phi~,al & al~,chi & chi~; 

6 input: TETRAD comprises 

• • • input: component I 0=E**al d phi/ sqrt (2) , 
8 input: component 1 1=C . C . (component I 0) , 

• • • input: component I 2=d x/ sqrt (2) , 

Winput: component I 3= (d y +2 RE(chi d phi) +h d x)/sqrt(2); 
11 input: end of data. 

Here the symbol 'al' stands for a (see Eq. (22)), 'C.C is the symbol of the complex 
conjugating operator while the origin of the other symbols seems to be evident. Continuing 
with a matter of notations, it is also useful to mention that here and below the enumeration 
at left hand side in the copies of fragments of GRGec input scripts and output listings, as 
well as the marks at left 'input: 1 (line of input code), 'output: 1 (line of output file), etc., 
are not a part of the texts displayed. They were introduced for convenience of further 
references alone. For example, we can refer to the lines 1,3,5,6-10,11 involved in the DATA 
section displayed above.) 

It is seen that the above script encodes, almost verbatim, the content of Eqs. (21, 22, 
23). Nevertheless a number of comments, concerning some specific points of GRGec input 
language, should be made. 

First, the tilde symbol '"' attached at right to some identifiers imitates the overscoring 
used to denote complex conjugation in the standard system of mathematical notations. 
For example, 'phi~' may be regarded as the representation of 0, etc. However, contrary 
to the habit of mathematical notations, the additional tilde mark does not impose itself 
any actual relation between the corresponding identifiers. For example, phi and phi" are 
a priori not connected in any way. The encoding of the necessary relationship such as 
C.C. ((f)) = is realized by means of a separate declaration. See line 5. 

Next, the term 'SCALARS' 13 denotes, roughly, a sort of functions each of whose is at the 
same time a single symbol (identifier). The dependence upon the (fixed) set of arguments 14 
pointed out in the declaration (line 3) is supported by a special implicit mechanism. 

12 One may organize a library of such sort data and load them from there. 
13 It has no relation to the own Reduce' scalars. 

14 which, to be more precise, is varied by the coordinate transformations, if any 
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We also mention that the symbol df (or, equivalently, DF) denotes partial derivative 
(as well as ordinary derivative of a function of a single argument which is not formally 
distinguished from a partial derivative). In particular, df (chi , phi) encodes dx/d(f). The 
derivative d 2 x/d<p/d(f> is represented by the script df (chi ,phi ,phi~), for d 2 x/d(p 2 one 
has to write down df (chi ,phi , 2) , etc. (These notations are identical to ones supported 
by Reduce system.) 

Finally, the record ' component I (j)' (the delimiter I is here optional and might be 
dropped out), where (j) is a digit, refers to a single component of a data object, in our 
case TETRAD. Thus the record 'component! 1= C.C. (component 1 0)' (see line 8) represents 
the equation 9° = ¥, cf. Eqs. (3), (4). 

Now we are ready to proceed with calculations. 

Let GRGec package was started and has read the problem specification consisting 
of a single section of DATA. We may to immediately control the calculations issuing the 
instructions from a keyboard. (Alternatively, they — or a part of them — might be put 
into a section of INSTRUCTIONS, provided it were included in the problem specification. 
The mixed form of control is also available.) 

Thus, at first, after the auxiliary request to 
input: turn on displaying of negative powers 15 
affecting the format of output, we issue the instruction to 
input: find UNDOTTED CONNECTION 

The system response to it comes to the following output: 
output: ==> find UNDOTTED CONNECTION 

output: UNDOTTED CONNECTION is not found in the DATA section... 

output: ...UNDOTTED CONNECTION has been calculated 1.0 s. 

output: Total time spent amounts to 3.7 seconds. 

(Remark: Hereinafter, the symbol denotes some blank space cut off from the copy of a 

line of an input script or an output listing in order to fit it to the page width.) 

However, we want to look at the result (UNDOTTED CONNECTION) in explicit form. To 
that end, we have to initiate the action 'TYPE' directed to the data object we are interested 
in. The shortened form of the corresponding instruction looks as follows: 
input: UGAMMAO? 

Here UGAMMAO is the reference to the component T of the UNDOTTED CONNECTION ('UGAMMA' 
is its kernel name while '0' is the value of the summed undotted index distinguishing the 
component, see section 2.1). The return of the instruction reads: 

output: ==> UGAMMAO? 

output: — > TYPE UGAMMAO 



output: 1 - al~ 

output: UGAMMA = ( *E *DF(chi~,y)) d x + ( 

output: 2 

output: 1 al 

output: - ( *E )*(DF(al,y) + DF(al~,y))) d phi 

output: 2 



which may be shortened to 'turn on DIV 
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Thus Tq does not automatically vanish. This simply means that the expression for the 
TETRAD assumed is more general than we actually need. At the same time the condition 
ensuring the vanishing of Tq is very simple: al+al~ (i.e. a + a) and chi~ (i.e. x) have 
to be independent on y (y). Additionally, we have seen above that the value of al-al~ is 
completely gauge dependent. In particular, it always may be assumed to be independent 
on y and, then, both al and al~ turns out to be independent on y as well. Besides, x is 
also independent on y. We 'inform' the system of these facts by means of the following 
instruction: 

input: Let DF (al , y) =0 , DF ( al~ , y ) =0 , DF ( chi , y ) =0 , DF ( chi ~ , y ) =0 

The response reads 

output: ==> Let DF(al, y) = 0, DF(al~, y) = 0, DF(chi, y) = 0, DF(chi~, y) = 

output: Section of SUBSTITUTIONS is not found. 
output: It is added to the problem specification. 

output: The following substitution rules 

output: (1) DF(al , y) -> 

output: (2) DF(al~, y) -> 

output: (3) DF(chi, y) -> 

output: (4) DF(chi~, y) -> 

output: are put to the section of SUBSTITUTIONS. 

output: Substitutions 1,2,3,4 are now active 
Now we can inspect the resulting structure of UND0TTED CONNECTION. The instruction 
38 input: evaluate UGAMMA0 , UGAMMA 1 / \d x/\d phi/\d phi~ ,UGAMMA2/\d x/\d phi- 
yields in particular the output 
output: The expression 
output: UGAMMA0 

output: vanishes. 0.7 s. 

output: The expression 

43 output: UGAMMAlA d x/\ d phi/\ d phi~ 

44 output: vanishes. 0.8 s. 

output: The expression 

46 output: UGAMMA2A d x/\ d phi~ 

47 output: vanishes. 1*^1 0.8 s. 
output: Total time spent amounts to 4.5 seconds. 

Thus the condition distinguishing radiative electrovacs (Eq. (15)) is now satisfied. Addi- 
tionally, we see that the holonomic expansion of UGAMMA 1 does not involve d y-term while 
the expansion of UGAMMA2 consists of the d x- and d phi~-terms alone. We shall make use 
of these facts below. (We might display the explicit expressions of all the components of 
UND0TTED CONNECTION but these would be, essentially, useless here). 

Our next operation implies the inferring of some information from the field equations. 
Specifically, let us ask the system to 
49 input: obtain and isolate 

■■■input: the SCALAR PART of EINSTEIN-MAXWELL EQUATIONS; 
51 input: show EQUATIONS 

The answer is the following: 

output: ==> obtain and isolate the SCALAR PART of EINSTEIN - MAXWELL EQUATIONS 

output: — > OBTAIN SCALAR PART of EINSTEIN - MAXWELL EQUATIONS 
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output: . . . SCALAR PART OF EINSTEIN - MAXWELL EQUATIONS < 

>has been obtained 4.0 s. 

output: Total time spent amounts to 8.6 seconds. 

output: — > ISOLATE SCALAR PART of EINSTEIN - MAXWELL EQUATIONS 

output: al al ~ 

output: Warning: The denominator E *E has been removed assuming it to be 
output: nonzero . . . 

output: EQ. number 1 is put to the list of EQUATIONS. 
output: Total time spent amounts to 8.7 seconds. 
output: ==> show EQUATIONS 

output: al + al~ 

output: Eq. (1) reads 0=4*(E *DF(h,y,2) - DF(al ,phi ,phi~) 

output: - DF(al~ ,phi ,phi~) ) 

Evidently, it makes sense to 

input: resolve EQUATION (1) w.r.t. df(h,y,2); 

The result is Eq. (2) (we shall display it explicitly below) while the preceding EQUATION 
may be removed (by means of the instruction 'erase EQUATION (1)'). 

The last calculation we carry out in frames of the current step is the deriving of a conse- 
quence of Eq. (2) (which determines df (h,y,2), see below) by means of its differentiating 
with respect to the COORDINATE y. Specifically, let us claim to 

input: isolate EQUATION DR(LHS(EQ(2))-RHS(EQ(2) ) ,y)=0; 

input: show EQUATIONS; 

{Remark: Here the symbol 'DR' is so called cross derivative operator which is similar the to 
partial derivative but takes into account the mutual tree- like dependences of the SCALARS, if any, 
regarding them as functions, not mere variables. Mathematically, one may regard DR(. . . ,y) in 
the above instruction as a partial derivative with respect to y, provided a proper interpretation 
of dependences of variables (SCALARS and COORDINATES) is implied.) This calculation yield a 
new EQUATION and we obtain in particular the following output: 

output: Eq. (3) reads DF(h,y,3)=0 

output: Eq. (2) reads DF(h,y,2) 

output: - al - al~ 

output: =E *(DF(al,phi,phi~) + DF(al~ ,phi ,phi~) ) 

The current point is just reasonable for the terminating of the present (the first) step 
of the problem solving. The result having been derived is the following: we determined 
the dependence of all the unknowns on the COORDINATE y (which represents the affine 
parameter y along the multiple null principal congruence). Specifically, the SCALARS al, 
chi are independent on it (as well as al~, chi~) while h is a second order polynomial in 
y, the higher order coefficient being found in explicit form. 

Additionally, we have seen that a gauge transformation allows to add any function of 
x, <p, to Xi4>+Xi<t>- See Eq. (26). Thus since both x an d X are independent on y the 
above 2-term sum can be nullified by means of a proper gauge choice. Then, there is a 
real 'potential' k = k(x,(fi,(j)) such that x — — ~~ i Ultimately, noting that the 

functions XiX enter the tetrad in the form of the expression 3?(xd0) alone, we have the 
reduced representation 

»(xd0 = -9(fc*d$ (27) 
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which will be used below in 'formulae' (GRGec scripts) determining TETRAD. 

An additional remark is now in order. It has to be mentioned that in this section 
we exhibited, essentially, the complete session of the work with GRGec- Indeed, the 
input code used in the first step of calculations consists, at first, of the section of DATA 
which was displayed in the beginning and, at second, a series of instructions (assumed 
here to be issued from a keyboard) which was listed (and commented) in the course of 
the discussion as well 16 . We also reproduced all the substantial output which follow the 
instructions execution (the only omission is the initialization messages and the outcome 
of some preliminary analysis of the initial data by the system including the copy of the 
structured problem specification). 

Unfortunately, in the sequel, in many cases the authentic copies of output would occupy 
too much space to be given here. For example, the output of the only next (the second) 
step of calculation considered below is more that 30 kB long, i.e., would require at least 10 
pages even in a condensed format. Hence we shall carry out below a strict selecting and 
shall display only those output which is actually necessary. On the contrary, the input will 
be displayed without omissions. 

3.2 Step 2: Introduction of electromagnetic field and further 
reduction of field equations 

The DATA incorporated in the problem specification of the second step of calculations 
develop ones used in the preceding step (see above lines l-ll). In particular, the symbols 
of COORDINATES need not be varied, of course. Similarly, the TETRAD is described by the 
same paragraph (a fragment of a section bounded at right by the semicolon ';') which 
comprises lines 6-10), where however the equation (27) is taken into account. Thus now 
instead of the line 10, one has the following one: 

73input: component I 3= (d y -2 IM(df (k,phi) d phi) +h d x)/sqrt(2); 

Furthermore, now the identifier h is not an abstract 'SCALAR'. Instead, it is specified an 
ABBREVIATION by means of the statement 

74 input: ABBREVIATION follows: h=R+Q*y+P*y**2 ; 

(Remark: An identifier of ABBREVIATION is immediately replaced by the corresponding expression 
associated with it whenever it is met in a mathematical formula). Thus we realize the form of 
the dependence of h upon y (second order polynomial) inferred in the preceding section 
from the field equations. The thee symbols R,P,Q are the new SCALARS instead of the single 
h but, as opposed to h, they are independent on y. We shall need even more SCALARS. 
Their total collection looks now as follows: 

75input: declare SCALARS al,k,R,P,Q(x,phi,phi~) , 

••• input: gamma, tau,zeta, lam,mu(x,y ,phi ,phi~) , 

77 input: psi , ems (x , y , phi , phi") ; 

(c/. line 3). Besides, let us 

16 To be completely precise, the file slang containing the additional section of CONFIDENTIAL SLANG which 
installs the extended collection of 'keywords' is read and included in the problem specification (cf. the list 
line of the script displayed in Appendix A). We omit discussion of this facility here however. See [1]. 
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78 input: declare REAL x,y,k,h,R,P,Q; 

input: Declare COMPLEX CONJUGATED phi & phi~,al & al~,lam & lam", 
input: mu & mu~,zeta & zeta~, 

input: psi & psi", ems & ems"; 

Further, we may specify the structure of UNDOTTED CONNECTION basing on its properties 
derived in the preceding subsection. Specifically, 
UNDOTTED CONNECTION comprises 

83 input: component I l=gamma d phi + tau d phi" + zeta d x, 

84 input: component I 2= (mu d phi" + lam d x)/E**al; 

(Remark: Dropping out specification of the component|0 of UNDOTTED CONNECTION (Tg in the 
mathematical notation) in this working binding of the data object value, it is automatically 
nullified.) This representation of connection components merely expresses the form of their 
expansions implied by the equations displayed in the lines 43-44, 46-47, see the comment 
following them. The denominator E**al (i.e. e a ) of the component|2 is introduced for the 
sake of the later convenience (it makes some relationships simpler). The SCALARS gamma, 
tau, zeta, mu, lam are assumed to be subjected no a priori constraints, being therefore 
defined just by the expansions displayed in lines 83-84. 

Accomplishing the DATA section of the current version of the problem specification, we 
add the following description of the electromagnetic field: 
input: UNDOTTED EM SPINOR components are 

86 input: component I l=psi , 

87 input: component I 2=ems ; 

88 input: DOTTED EM SPINOR is HERMITEAN CONJUGATED to UNDOTTED EM SPINOR; 

Since the SCALARS psi, ems still correspond to arbitrary functions (see line 77), it, essen- 
tially, only states that the componentlO (encoding <pg = <p 00 ) vanishes, implying addition- 
ally (by virtue of the statement displayed line 88) that </>^ vanishes as well. It is worth 
noting that the vanishing of <f)g (4> g ) is neither a restriction nor an ansatz since it can be 
deduced from Eqs. (19), (15) and the spinor part of Einstein-Maxwell equations (13) in a 
way similar to one used for the derivation of Eqs. (20). 

Proceeding now with calculations, let us first of all examine implications of the first 
structural equations. Since, given the tetrad, they have to uniquely determine the connec- 
tion forms, it may be expected that, in the case under consideration, the explicit values of 
the SCALARS gamma, tau, zeta, mu, lam have to follow. This is indeed the case. Namely, 
the instruction 
input: obtain and isolate 

input: the SECOND VERSION of the UNDOTTED FIRST CARTAN EQUATIONS 

entails just the five EQUATIONS (automatically enumerated with the numbers 1,. ..,5). 
Then, requiring to 

91 input: resolve EQUATIONS (l)-(5) w.r.t. gamma, tau, mu, zeta, lam 

(after that one may 'erase EQUATIONS (l)-(5)') and then, issuing the instructions 

92 input: turn on the displaying of negative powers; 

93 input: show new EQUATIONS factoring out E**al,E**al~ 

we obtain the following output: 

output: 1 

output: Eq. (10) reads gamma= *DF(al~,phi) 

output: 2 
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output: 1 

output: Eq. (9) reads tau= *DF(al,phi~) 

output: 2 

output: 1 al + al~ 

output: Eq. (8) reads mu= - ( *E )*(DF(al,x) + DF(al~,x)) 

output: 2 

output: - DF(k,phi,phi~)*I 

104 output: 11 1 

■■■output: Eq. (7) reads zeta= *DF(al,x) *DF(al~,x) + P*y + *Q 

• • • output: 4 4 2 

■■■output: 1 al~ -1 al -1 

■■■output: + *(E ) *(E ) *DF(k,phi,phi~)*I 

109 output: 2 

output: Eq. (6) reads lam= - DF(k,x,phi)*I - 2*DF(k,phi) *I*P*y 

output: 2 
output: - DF(k,phi)*I*Q + DF(R,phi) + DF(P,phi)*y 

output: + DF(Q,phi)*y 



These relationships are just the content of the first structural equations. Essentially, they 
define the connection forms, see lines 83-84. 

We may regard the expressions determining auxiliary SCALARS gamma, tau as the ulti- 
mate formulae and use them for the complete eliminating these SCALARS (which are involved 
in the expansions of UNDOTTED CONNECTION alone). To that end, we have preliminarily to 

input: add substitution rules LHS(EQ(9))->RHS(EQ(9)) , 

input: LHS(EQ(10))->RHS(EQ(10)) 

to the section of SUBSTITUTIONS. (Although we have not included this section it in the 
problem specification, it is automatically introduced.) (Remark: The above record of sub- 
stitution rules mean that, bringing them into action, the expression at left in EQ[.](9) displayed 
above is replaced by the expression situated at right (which is indicated by the symbol '->'); the 
second formula involving EQ[.](10) possesses similar meaning.) 

The substitutions introduced are automatically endowed with the numbers 1 and 2, 
respectively 17 . They are used for the referring to the rules. It particular, we invoke the 
substitutions which ensure the eliminating of gamma and tau in favor of their values by 
means of the instruction 

input: match the rules (1) , (2) with UNDOTTED CONNECTION 

(In order to verify that the necessary transformation has been indeed carried out one may 
issue the instruction 'type UNDOTTED CONNECTION'.) 

Next we consider a component of Einstein field equations which was also used in the 
first step of calculation. The instruction 

input: obtain and isolate the SCALAR PART of EINSTEIN-MAXWELL EQUATIONS 
entails the only new EQUATION endowed with the number 11. Requiring to 'show new 
EQUATION', we obtain in particular the output 
output: al + al~ 

output: Eq. (11) reads 0=2*E *DF(zeta,y) - DF(al ,phi ,phi~) 

output: - DF(al~ ,phi ,phi~) 

17 One might introduce any other numerical labels instead. 
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Notice that this EQUATION involves the derivative of the SCALAR zeta whose value has been 
specified by the above Eq. (7). A natural consequence of these two EQUATIONS, implying 
the elimination of DF(zeta,y), may be derived by means of the instruction 
input: isolate EQUATION LHS(EQ(11) )=RHS(EQ(11) ) 

input: -COEFFN (RHS (EQ ( 11) ) , df (zeta , y) ) *DR (LHS (EQ (7) ) -RHS (EQ (7) ) , y) 
This is an example of operation frequently used for reductions of overdetermined systems 
of equations with partial derivatives. Removing the second line we would obtain the EQUA4- 
TION identical to the initial one (number 11). The additional term in the lower line is 
proportional to the derivative of the discrepancy 'LHS(EQ(7))-RHS(EQ(7))' of Eq. (7). 
Thus the resulting EQUATION (number 12) is equivalent to Eq. (11) modulo Eq. (7). The 
coefficient in front of the derivative of the discrepancy is chosen in such a way to automati- 
cally entail the desirable simplification. In our case we have to eliminate df (zeta,y); hence 
it equals COEFFN (RHS (EQ( 11) ) ,df (zeta.y) ), i.e. the coefficient in front of the derivative 
df (zeta,y) involved in r.h.s. expression of Eq. (11) which is considered here as polynomial 
in df (zeta,y) (a linear function in fact). 

It turns out that the last EQUATION obtained determines the value of the SCALAR P. The 
instructions 

input: resolve Eq. (12) w.r.t. P; show new EQUATIONS 
yield in particular 
124 output: 1 - al - al~ 
■■■output: Eq. (13) reads P= *E 

• • • output: 2 

127 output: *(DF(al,phi,phi~) + DF(al~ ,phi ,phi~) ) 

which will play an important role in what follows. 

Now the optimal way of the further reduction of the field equation is to proceed with 
Maxwell equations. Thus the next instruction reads: 
input: obtain and isolate UNDOTTED MAXWELL EQUATIONS 

It yields four 'plain' EQUATIONS (number 14-17) which can be found to determine a number 
of SCALAR derivatives. Specifically, the instructions 

129 input: resolve Eqs . (14) -(17) w.r.t. 

130 input: DF(psi ,y) ,DF(ems ,y) ,DF(psi ,phi~) ,DF(psi ,phi) ; 
input: show new EQUATIONS 

yield the corresponding explicit expressions (Eqs. (18)-(21)). We obtain in particular 
that 

132 output: Eq. (21) reads DF(psi,y)=0 

and 

133 output: Eq. (19) reads DF(psi ,phi~)=0 

(We do not display the values of DF(ems,y) (Eq. (20)) and DF(psi,phi) (Eq. (18)) since 
it is only essential currently to know that these are nonzero). We see therefore that 
the function if) (represented by the SCALAR psi) which equals the component 4>i of the 
undotted spinor of electromagnetic field (see line 86) depends on the coordinates x and 
alone (represented by the symbols x, phi). Since is a complex variable ip has to be 
holomorphic with respect to it (and is assumed to be C°° with respect to x). Further, 
we have to inform the algebraic processor that the constraints 132,133 are to be taken 
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into account throughout all the subsequent calculations. This is realized by means of the 
instruction 

input: let the last EQUATION hold true, 
input: the last but 2 EQUATION hold true 

A series of the further manipulations with the collection of EQUATIONS obtained so far 
follows a standard routine procedure which widely applies for the reducing of overdeter- 
mined systems of quasilinear equations with partial derivatives 18 . Specifically, some of the 
equations available are differentiated and, then, combined in such a way to cancel out the 
higher order derivatives arisen. Such an operation is repeated until it yields new tractable 
relationships. Thus, following this scheme, the concrete transformations, which lead to the 
desirable simplifications and entail new data (new EQUATIONS), are chosen by ourselves 
with the help of the analysis of the intermediate output. In all the cases we meet the latter 
is fairly straightforward but, nevertheless, too bulky to be displayed here. Hence we shall 
not exhibit all these speculations and shall outline only the most substantial points of the 
derivation whose code is displayed below. 

Specifically, the following series of instructions is issued: 

136 input: isolate EQUATION DR(LHS(EQ(18) ) ,y)=DR(RHS(EQ(18) ) ,y) ; 

137 input: show new EQUATIONS; 

138 input: resolve Eqs . (20) , (23) w.r.t. DF(ems,y) ,DF(ems,y,phi~) ; 

139 input: show new EQUATION; 

UOinputlet (10) the last EQUATION hold true; 

141 input: isolate EQUATION DR(LHS (EQ (20) ) ,y)=DR(RHS(EQ(20) ) ,y) ; 

142 input: resolve Eqs. (25) , (20) w.r.t. DF(ems ,y , 2) , ems ; 

143 input: abolish substitution (10); 
input: show Eq. (27) ; 

145 input: let the last EQUATION hold true; 

146 input: resolve Eqs. (18), (26) w.r.t. DF(k,phi,phi~) ,df (ems,y) ; 
input: show new EQUATIONS; 

148 input: resolve Eqs. (26) , (28) w.r.t. df (ems ,y) , ems ; 
input: show new EQUATION; 

Eq. (18) (produced by the instruction displayed in lines 129-130) expresses DF(psi,phi). 
By virtue of Eq. (21) L.H.S. of the EQUATION described in the line 136 vanishes. Explicitly, 
the corresponding new 

input: Eq. (23) reads 0=DF(al,phi~)*DF(ems,y) + DF(ems ,y ,phi~) 
(a part of the output of line 137). In turn, Eq. (20) expresses DF(ems,y) and lines 138- 
139 yield separate representations for DF(ems,y,phi~) and DF(ems,y). The instruction 
yields a single new EQUATION (expressing DF(ems,y,phi~) 19 ). Line 140 initiates subsequent 
replacing of the derivative DF(ems,y,phi~) by the corresponding value, the rule being 
endowed with number 10. (In line 143 this substitution is disabled.) Eq. (20) determines 
DF(ems,y). Thus in the line 141 the second order derivative DF(ems,y,2) is introduced. 

The instruction in line 142 eliminates ems from R.H.S. of the new EQUATION (number 
25). One of the resulting EQUATIONS is most important, namely, 
151 output: Eq. (27) reads DF(ems,y,2)=0 

18 For the first reading it might be recommended to skip the fragment below up to the line 151. 
^Representation of DF(ems,y) does not vary; its inclusion to the instruction 138-139 is intended for the 
eliminating of DF(ems,y). Such a trick will be often used below. 
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Thus ems is linear in y. Line 145 ensures the nullifying of DF (ems , y , 2) in all the subsequent 
calculations. 

It is straightforward now to determine the coefficient in front of y in ems expansion. 
The instruction in line 146 yields a pair of EQUATIONS. In particular 

152 output: - al 

153 output: Eq. (28) reads DF(ems,y)=E *DF (psi ,phi) 

(psi does not depends on y). Finally, line 148 yields representation of ems through 
DF(ems,phi~) and the expressions independent on y (Eq. (30)) which will be utilized 
later on. 

Next, we shall deduce some useful consequences of Einstein-Maxwell equations. To 
that end, let us, at first, 

input: obtain the SPINOR PART of EINSTEIN-MAXWELL EQUATIONS; 
input: isolate the component 1 1 1 1" of the ABOVE EQUATIONS; 
input show new EQUATIONS 

This yields in particular the output: 
output: al + al~ al + al~ 

output: Eq. (31) reads 0= - 2*E *DF(zeta,y) - E *psi*psi~ 

output: - DF(al,phi,phi~) - DF(al~ ,phi ,phi~) 

Eq. (31) involves the derivative (zeta,y). On the other hand we had the 'explicit' represen- 
tation of zeta (Eq. (7), see lines 104-109). The dependence of zeta on y is characterized by 
the SCALAR P which, in turn, is determined by Eq. (13) (lines 124-127). These relationships 
(determining DF(zeta,y)) can be represented by means of the two substitution rules which 
are introduced as follows: 

input: add substitution rules 

input: (30) P -> P-(LHS(EQ(13))-RHS(EQ(13))) , 

input: (31) DF(zeta,y) -> DF (zeta, y) -MATCHING (DR(LHS(EQ (7) )-RHS(EQ (7)) ,y) ,30) 
Notice that, formally, the rule (30) is trivial modulo Eq. (13). However, the term involving 
it is chosen in such a way to eliminate P at right. Similarly, the rule (31) is trivial modulo 
some expression vanishing if Eq. (7) is satisfied. {Remark: Here the macro MATCHING applies 
the substitution rule whose number is given as the second argument 20 (i.e., the rule (30)) to the 
first argument.) The combined effect is the eliminating of the both DF(zeta,y) and P at 
right in the expression to be substituted in accordance with the rule (31). Accordingly, 
applying it to Eq. (31) by means of the instruction 

input: isolate EQUATION 0=MATCHING(LHS(EQ(31) )-RHS(EQ(31) ) ,31) ; 
(here '31' in EQ(31) is the number of an EQUATION while the last '31' is the number of a 
substitution rule; their coincidence is occasional), the request 'show new EQUATION' yields 
in particular the following output: 
164 output: al + al" 

■■■output: Eq. (32) reads 0=E *psi*psi~ + 2*DF(al,phi,phi~) 

166 output: + 2*DF(al~,phi,phi~) 

This EQUATION is worth a separate attention. It involves the SCALARS al, al", psi, 
psi~ (representing a,a,ifj,ijj) which may be regarded as functions of some fixed sets of 
arguments. In accordance with declaration shown in lines 75-77, al and al~ depend on x, 

20 An arbitrary number the rules is allowed to be invoked by the 'MATCHING' operator. 
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phi, and phi~ while more strong restriction is imposed on psi (see lines 132, 133): it may 
depend on x and phi alone (i.e., is holomorphic with respect to phi). Then, since psi and 
psi~ are COMPLEX CONJUGATED, the SCALAR psi~ depends on x and phi~ (and represents 
the function of x and holomorphic with respect to the second argument). Taking these 
relationships into account, we see that Eq. (28) may be regarded as Lioville equation with 
respect to unknown al+al~ (2 Re a) for some given psi, psi~ (tp and ■0), holomorphic and 
antiholomorphic, respectively (cf. [9], Eq. (27.50)). Besides, we have seen above that the 
value of al-al~ is a matter of a gauge, provided it does not depend on the COORDINATE y. 
Other relevant gauge freedom involves the transformation of the COORDINATE phi (and, si- 
multaneously, its complex conjugated phi~) which, in mathematical notations, is described 
by the formula — > H(0, x) for arbitrary S holomorphic with respect to 0, assuming the 
restriction S(0, x),^^ to be fulfilled. The latter fact has to be taken into account when 
analyzing solutions of Lioville equation and their gauge transformations. 

It should be noted now that GRGec is not able (and was not intended) to solve differen- 
tial equations 21 . We have to perform this work ourselves. In our case it is straightforward 
to find and write down the general solution to Eq. (32) (cf. [9], Eq. (27.51)). It is conve- 
nient to represent the solution in the form of substitution rules. They are introduced by 
the following instruction: 

input: add rules 
I68input: (40) al -> -log(psi) -log(l+phi*phi~/4) , 
■■■input: (41) al~ -> -log(psi~)-log(l+phi*phi~/4) , 

■■■input: (42) df (al ,phi ,phi~) -> DR(-log(psi)-log(l+phi*phi~/4) ,phi,phi~) , 
nimput: (43) df(al~,phi,phi~) -> DR(-log(psi~)-log(l+phi*phi~/4) ,phi ,phi~) ; 

(here, evidently, the third and the fourth rules are the consequences of the first and second 

ones) which just describe the general local solution to Eq. (32)). 

It is instructive to show that we indeed deal with a solution the Lioville equation (Eq. 

(32) ). Preliminarily, we have to inform the system on the properties of the SCALAR psi~ (as 

opposed to psi, we have not dealt with it yet in fact). It is convenient to do this referring 

to the corresponding EQUATIONS concerning its complex conjugated counterpart, psi, i.e. 

Eq. (19) and Eq. (21). Taking their complex conjugation we obtain the appropriate 

characteristics of psi". To that end, let us issue the instructions 
input: evaluate aux=C . C . (LHS (EQ ( 19) ) ) ; 
input: let (50) aux->C . C . (RHS (EQ (19) ) ) ; 
mput: evaluate aux=C.C. (LHS(EQ(21))) ; 
input: let aux->C.C. (RHS (EQ (21) ) ) ; 

(Remark: Here and below the ('user-defined') identifier aux plays role of auxiliary variable used 
for the temporary storing the expressions situated on r.h.s. of assignments in instructions with 
the action EVALUATE.) These substitution rules are labeled with numerical labels 50 (specified 
by ourselves) and 51 (the next number generated automatically). 
Then the instruction 

176 input: evaluate aux=MATCHING (LHS (EQ (32) ) -RHS (EQ (32) ) , 42 , 43) , 

177 input: MATCHING (aux , 40 , 41) ; 

21 This point should not be considered as its inalienable feature. Moreover, it could be worth introducing 
an interface with an appropriate ODE or PDE package. However such a facility has not been realized so 
far. 
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reports zero value which just means the expectable satisfaction of the Eq. (32) (see lines 
164-166) by virtue of the relationships displayed in lines 168-171. (Let us remind that lines 
176,177 re-evaluate the 'discrepancy' of Eq. (32) executing, at first, substitutions number 
42,43, and then, substitutions number 40,41 which were specified in lines 168-171). 

An additional useful result immediately follows from Eq. (13) and Eq. (32). Their 
appropriate superposition is realized by means of instructions 



input: isolate EQUATION 

input: LHS(EQ(9))=RHS(EQ(9)) 

input: -C0EFFN (NUMR (RHS (EQ (9) ) ) , DF (al , phi , phi" ) ) 

input: /DENM(RHS(EQ(9)))/2 

input: * (LHS (EQ (280) ) -RHS (EQ (28) ) ) ; 

183 input: show equation (29); 

input: isolate EQUATION 

input: LHS(EQ(13))=RHS(EQ(13)) 

input: + (COEFFN (NUMR (RHS (EQ ( 13) ) ) , DF (al , phi , phi~ ) ) 

input: +C0EFFN (NUMR (RHS (EQ ( 13) ) ) , DF (al~ , phi , phi" ) ) ) 

input: /DENM(RHS(EQ(13) ) ) 

input: * (LHS (EQ (32) ) -RHS (EQ (32) ) ) 

input: I (COEFFN (NUMR (RHS (EQ (32) ) ) , DF (al , phi , phi" ) ) 

input: +C0EFFN(NUMR(RHS(EQ(32))) ,DF(al~ ,phi ,phi~) ) ) ; 

input: show new EQUATION; 

In particular, the last line entails the the response 

193 output: 1 

■■■output: Eq. (33) reads P= *psi*psi~ 

195 output: 4 



The last action we shall carry out in frames of the current step is the explicit separating 
of the dependence of the SCALAR ems on the COORDINATE y. As we have seen it is a linear 
function of y (see line 15l). Moreover, the coefficient in front of y is specified by Eq. (28) 
(lines 152-153). Thus it retains to find the free term which can be explicitly represented as 
ems-y*df (ems , y) . Instead we shall calculate the derivative of (ems-y*df (ems , y) ) *E**al 
with respect to phi". The corresponding instructions are straightforward, executing several 
transformations which are identical ones modulo Eq. (28) and Eq. (30): 
input: evaluate 

197 input: expr= (ems-y*df (ems ,y) ) *E**al 

198 input: aux=DR(expr,phi~) , 

input: aux=aux-COEFFN(NUMR(aux) , df (ems , phi" ,y) ) /DENM(aux) 

input: *DR(LHS(EQ(28))-RHS(EQ(28)) ,phi~) , 

input: aux=aux-COEFFN(NUMR(aux) , df (ems ,y) ) /DENM(aux) 

input: * (LHS (EQ (28) ) -RHS (EQ (28) ) ) , 

input: aux=aux-C0EFFN (NUMR(aux) , ems) /DENM (aux) 

input: * (LHS (EQ (30) ) -RHS (EQ (30) ) ) , 

input: aux 

input: factoring out E**al , E**al~ ; 

Here the line 197 plays role of the definition of expr to be estimated, line 198 introduces its 
derivative denoting it aux which is further undergone transformations 'identical' modulo 
Eq. (28), Eq. (30). The last two lines of the code exhibits its resulting representation: 

output: The expression 

output: aux 

output: amounts to: 
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output: al + al~ 

output: - E * (DF(al ,x) *psi + DF(psi,x) + DF(al~ ,x) *psi) 

output: - 1 

output: + I *(2*DF(k,phi,phi~)*psi + DF(k,phi~) *DF(psi ,phi) ) 1<~~>1 1.3 s. 



The dependence derived can be partially integrated by means of the extracting from the 
(value of) auxiliary variable expr a series of appropriate expressions which, after the sub- 
tracting their corresponding derivatives from the value of aux (equal to DR ( expr , phi ~) 
modulo Eq. (28) and Eq. (30)), lead to the cancelling out certain terms of the latter. 
Specifically, the following instructions are to be issued: 



input: evaluate 

IXhinput: expr=expr -df(k,phi)* COEFFN(NUMR(aux) , df (k,phi , phi") ) /DENM(aux) , 

• • • input: aux=DR ( expr , phi ~ ) , 

■■■input: aux=aux-COEFFN(NUMR(aux) ,df (ems ,phi~ ,y) )/DEM(aux) 
-input: *DR(LHS(EQ(28))-RHS(EQ(28)) ,phi~) , 

■■■mput: aux=aux-COEFFN(NUMR(aux) ,df (ems ,y) ) /DEM (aux) 
■■■input: *(LHS(EQ(28))-RHS(EQ(28))) , 

■■■input: aux=aux-COEFFN(NUMR(aux) ,df (k, phi , phi") ) /DEM (aux) 
222 input: * (LHS (EQ (29) ) -RHS (EQ (29) ) ) , 

input: aux ; 

input: evaluate 

225 input: expr=expr -k* COEFFN(NUMR(aux) , df (k,phi~) ) /DENM(aux) , 

• • • input: aux=DR ( expr , phi ~ ) , 

■■■input: aux=aux-COEFFN(NUMR(aux) ,df (ems ,phi~ , y) ) /DEM (aux) 
■■■input: *DR(LHS(EQ(28))-RHS(EQ(28)) ,phi~) , 

■■■input: aux=aux-COEFFN(NUMR(aux) ,df (ems ,y) ) /DEM (aux) 
■■■input: *(LHS(EQ(28))-RHS(EQ(28))) , 

■■■input: aux=aux-COEFFN(NUMR(aux) ,df (k, phi , phi") ) /DEM (aux) 
232 input: * (LHS (EQ (29) ) -RHS (EQ (29) ) ) , 

input: aux ; 

input: evaluate expr , aux-DR(-E** (al+al~) *psi ,x) 
input: factoring out E**al; 



Here the first application of the action 'EVALUATE' (lines 215-222) 'integrates out' the term 
proportional to df (k,phi) while the second application (lines 225-232) makes the same 
thing with the term proportional to k. The output of the last instruction looks as follows: 



output: The expression 
output: expr 

output: amounts to: 

239 output: al 
■■■output: E *( - DF(ems,y)*y + ems) 

• • • output: - 1 

242 output: +1 *( - 2*DF(k,phi) *psi - DF(psi,phi)*k) 0.2 s. 

output: The expression 

output: aux - DR(- E** (al + al~) * psi, x) 

output: vanishes. 0.2 s. 



Thus the derivative of the expression displayed in lines 239-242 with respect to phi" equals 
the derivative of -E** (al+al~) *psi with respect to x. Hence we obtain the free (y- 
independent) term of ems by a quadrature. 

Now the amount of the new information deduced suffices to renew the problem speci- 
fication. 
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3.3 Step 3: Accomplishing of reduction of the set of field equa- 
tions 

We modify the preceding problem specification is such a way to implement all the rela- 
tionships obtained so far. We need now to 
input: Declare SCALARS k,R,Q(x,phi,phi~) , 



247 input: psi , beta (x, phi) ,psi~ (x,phi~) , 

248 input: xi(x,phi,phi~) , 
input: rho (phi , phi ~ ) , 

input: lam(y , Q ,psi ,psi~ , x,phi ,phi~) , 

input: zeta(y,Q,rho,psi,psi~,x,phi,phi~) , 

input: mu(rho,psi,psi~,x,phi,phi~) , 

input: emsO(rho,psi,k,xi,x,phi,phi~) , 

input: emsl(rho,psi,x,phi) ; 



In addition to the SCALARS used before and displayed in lines 75-77 some new ones are 
introduced. The collections of their 'arguments' are determined by the VALUES which has 
to be bound with them. The latter realize the relationships (equations) derived above. 
Specifically, 

input: SCALAR VALUES follow: 



256 input: rho=l +phi*phi~/4, 

input: lam=DR(P,phi)*y**2 +y* (df (Q , phi) -2i*P*df (k,phi) ) 
input: -i*df (k,x,phi) -i*Q*df (k.phi) +df(R,phi), 

input: zeta=P*y +DR(al-al~ ,x)/4 

input: +0/2 +df (k,phi,phi~)*i/2/E**(al+al~) , 

input: mu= -DR(E**(al+al~) ,x)/2 -i*df (k,phi ,phi~) , 

input: emsl=E** (-al) *df (psi , phi) , 

input: emsO=E** (-al) * (-2i*psi*df (k,phi) -i*k*df (psi ,phi) 
input: -xi -df (beta, phi) ) ; 



As opposed to the interpretation followed to in the preceding steps, the symbols al (to- 
gether with al~) and P are now subsidiary objects. They are defined in accordance with 
the following statement: 

265 input: ABBREVIATIONS comprise 

••• input: al=-log(rho) -log(psi) , 

■■■input: P=-psi*psi~/4, 
268 input: h=R+Q*y+P*y**2 ; 

Cf. lines 168, 193-195, 256. (Here h definition exactly replicates the one used above, see line 

74). 

The declaration of REAL objects (see line 78) is now to be replenished with the SCALAR 
rho while all the other new SCALARS are defined as complex (by means of the declara- 
tion 'COMPLEX CONJUGATED'), being endowed with complex conjugated counterparts whose 
identifiers are marked with the character '"' attached at right. 

The specification of the basic geometric geometric objects — TETRAD and UNDOTTED 
CONNECTION — consists of the same formulae displayed in lines 6-10 and 83-84. As to the 
UNDOTTED EM SPINOR, its component 1 1 looks as above while the component 1 2 is endowed 
with the expansion revealing its dependence on y: 
input: component | 2=ems0+y*emsl ; 

This line replaces line 87. In particular, the former SCALAR ems (depending on y) is now 
replaced by the two new ones, emsO and emsl, which are independent on y. 
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The last — and novel — element of the problem specification is the following section 



of 






input: 


Substitutions : 


input: 


(1) 


lam->VAL(lam) ; 


input: 


(2) 


zeta->VAL(zeta) ; 


input: 


(3) 


mu->VAL(mu) ; 


274 input: 


(4) 


rho->VAL(rho) ; 


input: 


(5) 


emsl->VAL(emsl) ; 


input: 


(6) 


emsO->VAL(emsO) ; 


277 input: 


(7) 


df (xi,phi~)->DR(e**(al+al~)*psi,x 


input: 


(8) 


emsl~->VAL(emsl~) ; 


input: 


(9) 


emsO~->VAL(emsO~) ; 


input: 


end 


of substitutions. 



These list of substitution rules will be utilized whenever necessary. All they but one 
(displayed in line 277) mean the replacing of the identifier of a SCALAR by the corresponding 
SCALAR VALUE ('extracted' by the macros VAL). The specific substitution situated in line 
277 characterizes xi as the integral of the expression at right hand side with respect to the 
variable (complex COORDINATE) phi". 
Now let us proceed with calculations. 

At first, it is instructive to demonstrate how one can check that the current initial data 
verifies the part of field equations considered during the preceding steps of our calculation. 
This is carried out by means of the following instructions: 

input: obtain the SECOND VERSION of the UNDOTTED FIRST CARTAN EQUATIONS, 
input: UNDOTTED MAXWELL EQUATIONS, 

input: and the SCALAR PART of EINSTEIN-MAXWELL EQUATIONS; 

284 input: match substitution rules (1) - (4) , (5) - (7) with 

285 input: the ABOVE EQUATIONS; 

286 input: renew and type 

input: the ABOVE EQUATIONS; 

These equations are satisfied not automatically but as a consequence of specific SCALAR 
VALUES given in the problem specification. Accordingly, their substituting realized by the 
instruction displayed in lines 284-285 is necessary. The output of the above action TYPE 
situated in line 286 says us that all the equations listed are satisfied. Thus to obtain the 
'complete' solution of the problem it retains to 

input: obtain the SPINOR PART of EINSTEIN-MAXWELL EQUATIONS 

and ensure its fulfillment. 

In this case, similarly to the above procedure, we also have to take into account the 
SCALAR VALUES given (evaluating derivatives of SCALARS, they are automatically taken into 
account). To that end, the following instructions are to be issued: 
289 input: match substitution rules (l)-(3) , (5)-(6) , (8)-(9) with 

■■■input: the ABOVE EQUATIONS; 

291 input: isolate the ABOVE EQUATIONS; 

The result of the line 291 is four differential EQUATIONS (marked with the ordinal numbers 
1,2,3,4). The first of them is in fact trivial, being immediately fulfilled, provided the 
SCALAR VALUE for rho (see line 256) is taken into account. Indeed, the instruction 
input: Evaluate MATCHING (LHS(EQ(D)-RHS(EQ(l)) ,4) ; 
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reports zero result. (The substitution number 4 which is referred to here causes the re- 
placing of the SCALAR rho by the corresponding VALUE, see lines 274,256.) On the contrary, 
a single Eq. (4) entails two more EQUATIONS (which prove to be non-trivial) since it means 
the vanishing of a linear function of y. Accordingly, the instruction 
input: isolate 

input: EQUATION 0=COEFFN(LHS (EQ (4) ) -RHS (EQ (4) ) , y , 1) , 

input: EQUATION OCOEFFN (LHS (EQ (4) ) -RHS (EQ (4) ) , y , 0) ; 

constructs the EQUATIONS (their ordinal numbers are 5 and 6) which just express the van- 
ishing of the first (the first line, the last parameter equals 1) and zero (the second line, 
the last parameter equals 0) order coefficients. The resulting collection of four EQUATIONS, 
which all the field equations have been reduced to, can be resolved with respect to deriva- 
tives of unknown SCALARS and displayed by means of the instructions 
input resolve EQUATIONS (2) , (3) , (5) , (6) w.r.t. 

input: df (Q ,phi~) ,df (Q ,phi) ,df (R,phi ,phi~) , df (Q ,phi ,phi~) ; 

input show new EQUATIONS; 

The output looks as follows: 
output: 3 2 

output: Eq. (10) reads DF(Q,phi~) = ( - 2*DF(k,phi,phi~,2)*I*psi*psi~ *rho - 

<gc the continuation dropped out consists of 13 non-empty lines ^> 
output: Eq. (9) reads DF(Q,phi) = ( 

output: l<r^>l 2 2 

output: 2*DF(k,phi,phi~)*DF(psi,phi)*I*psi *psi~*rho 

the continuation dropped out consists of 13 non-empty lines ^> 

iOAoutput: Eq. (8) reads DF(R,phi,phi~) = ( 

output: l<*~*l 4 4 4 

output: 1<~~~>1 4*DF(k,phi,phi~,2)*DF(k,phi)*psi *psi~ *rho 

the continuation dropped out consists of 59 non-empty lines ^> 

output: Eq. (7) reads DF(Q,phi,phi~) = ( - DF(psi ,x) *psi~ 

output: l<r^>l 2 

output: ?<~~>? + DF(psi,phi)*DF(beta~,phi~)*psi*psi~*rho 

the continuation dropped out consists of 9 non-empty lines ^> 

These EQUATIONS express the content of all the field equations which have not been satisfied 
so far. Their further reduction is now in order. 

Specifically, we have three equations determining derivatives of the SCALAR Q. Thus they 
have to satisfy some consistency conditions. Calculating 2*RHS (EQ (7) ) -DR(RHS (EQ (9) ) , 4. 
phi~)-DR(RHS(EQ(10) ) ,phi) (by means of the action EVALUATE), one obtains identical 
zero. This means that Eq. (7) does not yields a new information (and thus may be 
ERASEd). On the contrary, the condition of mutual consistency of Eq. (9) and Eq. (10) is 
nontrivial. We shall see it reduces in fact to the equation determining the SCALAR k. The 
optimal rout of its derivation is the following. 

Let us introduce the auxiliary variable 'expr' which initially is assigned with the value 
Q. Then we make a series of additive transformations of expr. Simultaneously we shall 
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calculate another auxiliary variable aux by the formula aux=DR( expr, phi), replacing here 
the derivative of Q by its expression provided by Eq. (9). The corresponding instructions 
read: 

input: evaluate expr=Q , 



input: expr=expr -df (k,phi ,phi~) 

input: ' *C0EFFN(NUMR(RHS(EQ(9) ) ) ,df (k,phi , 2 ,phi~) ) 

input: /DENM (RHS (EQ (9) ) ) , 

input: aux=DR(expr , phi) - (LHS (EQ (9) ) -RHS (EQ (9) ) ) , 

input: aux ; 

input: evaluate expr=expr -k*C0EFFN (NUMR (aux) , df (k, phi )) /DENM (aux) , 

input: aux=DR(expr , phi) - (LHS (EQ (9) ) -RHS (EQ (9) ) ) , 

input: aux ; 
input: evaluate expr=expr 

input: -df (ps i , x) *C0EFFN (NUMR ( aux) , df (ps i , x , phi ) ) /DENM ( aux) 

input: -beta*COEFFN (NUMR (aux) ,df (beta, phi) ) /DENM(aux) , 

322 input: aux=DR(expr , phi) - (LHS (EQ (9) ) -RHS (EQ (9) ) ) , 

input: aux ; 

The output of the last line reads: 

output: The expression 

output: aux 

output: amounts to: 

output: psi~*xi 

output: 0.6 s . 

output: 2 



In accordance with line 322 aux equals modulo Eq. (9) the derivative of expr with 
respect to phi. This dependence can be integrated. Let us notice also that expr possesses 
a complex value. The instructions 



input: turn on the displaying of negative powers, 

input: the support of complex numbers; 22 

input: evaluate (expr+C . C . (expr) ) /2 , (expr-C . C . (expr) ) /2/i ; 
yield the output: 

output: The expression 

output: (expr + C.C. (expr))/ 2 

output: amounts to: 

output: 1 -11 -1 1 

output: *DF(psi,x)*psi + *DF(psi~ ,x) *psi~ + Q *psi*beta~ 

output: 4 " ' 4 "' 4 

output: 1 

output: *beta*psi~ 0.2 s. 

output: 4 

output: The expression 

output: (expr - C.C. (expr))/ 2/ i 

output: amounts to: 

output: 2 1 -1 

output: - DF(k,phi,phi~)*psi*psi~*rho *DF(psi,x)*psi 

output: 4 



this one may be shortened to 'turn on DIV, COMPLEX' 
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output: I -11 I 

output: + *DF(psi~,x)*psi~ *k*psi*psi~ *psi*beta~ 

output: 4 2 4 

output: I 

output: + *beta*psi~ 0.2 s. 

output: 4 



Thus after the integrating mentioned above, the separation of the real and imaginary parts 
of the equation obtained obviously yields explicit representation of Q and DF(k,phi,phi~), 
respectively. 

This result accomplishes the current (third) step of our calculations. 



3.4 Step 4: Determination of R 

It was shown above that unknown Q and the derivative DF(k,phi,phi~) can be explicitly 
expressed in the form involving the integral of the SCALAR function xi (depending on 
x,phi,phi~) with respect to phi. It is reasonable therefore to introduce another SCALAR 
(we denote it int_xi) whose derivative with respect to phi equals xi. Another SCALAR 
independent of phi arises as a 'constant of integration'. We denote it gam", gam" may 
depend on x and phi~. (Thus it represents the function j(x,$) holomorphic with respect 
to <j>.) 

Accordingly, the lines 247 and 248 of the declaration specifying SCALAR dependences 
now has to look as follows: 



input: psi ,beta(x,phi) ,gam~ ,psi~ (x,phi~) , 

input: xi , int_xi (x,phi ,phi~) , 

while its other lines are not changed. Besides, the items 'gam & gam~, int_xi & int_xi~ 
are added to the description of 'COMPLEX CONJUGATED' objects which now reads: 
input Declare COMPLEX CONJUGATED phi & phi~,al & al~,psi & psi", 
input: lam & lam~,mu & mu~,zeta & zeta", 

input: beta & beta", gam & gam", 

input: xi & xi", int_xi & int_xi~, 

input: emsO & ems0~,emsl & emsl~; 

Further, the list of ABBREVIATIONS (see lines 265-268) is replenished with two new items: 

361 input: dd_k=-k/2/rho**2 
■■■input: +i/(4 rho**2*psi*psi~) 

• • • input: * (aux-C . C . (aux) , WHERE 

■■■input: aux=-df (psi ,x) /psi-psi* (beta~+int_xi~+gam) ) , 

■■■input: Q_val=(l/4)*(aux+C.C. (aux) , WHERE 
366 input: aux=-df (psi , x) /psi+psi* (beta~+int_xi~+gam) ) ; 



which represent the values of DF(k,phi,phi~) and Q, respectively, derived in the preceding 
subsection. 

We shall use below the following 
input: Substitutions : 



input: (1) lam->VAL(lam) ; 

input: (2) zeta->VAL (zeta) ; 

input: (22) zeta->VAL (zeta) ; 

input: (3) mu->VAL(mu); 

input: (33) mu~->VAL(mu~) ; 

input: (4) rho->VAL (rho) ; 
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input: 


(5) 


input: 


(55) 


input: 


(6) 


input: 


(66) 


input: 


(7) 


input: 


(77) 


input: 


(8) 


input: 


(9) 


input: 


(10) 


input: 


(11) 


input: 


(12) 


input: 


(13) 


input: 


(14) 


input: 


(15) 


input: 


(16) 


input: 


(17) 


input: 


(18) 


input: 


(19) 


input: 


(20) 


input: end of 



emsl->VAL(emsl) ; 

emsl~->VAL(emsl~) ; 

emsO->VAL(emsO) ; 

emsO~->VAL(emsO~) ; 

df (xi , phi~) ->DR(e** (al+al~) *psi , 

df (xi~ , phi) ->DR(e** (al+al~) *psi~ 

df (k,phi,phi~)->dd_k; 

df (k,phi,2,phi~)->DR(dd_k,phi) ; 



x); 
x); 



Q->Q_val; 

df (Q, phi , phi ~ ) ->DR(Q_val , phi , phi ~ ) ; 
df (Q,phi)->DR(Q_val,phi) ; 
df (Q,phi~)->DR(Q_val,phi~) ; 
df(int_xi, phi)->xi; 
df (int_xi~ ,phi~)->xi~ ; 
df (int_xi, phi ,phi~) ->df (xi ,phi~) ; 
df (int_xi~ ,phi ,phi~)->df (xi~ ,phi) ; 
xi->df (int_xi , phi); 
xi~->df (int_xi~ , phi~) ; 
substitutions . 

Their meaning is manifest. (Remark: Let us remind that the macro 'VAL' returns the SCALAR 
VALUE of its argument, a SCALAR identifier.) 

The aim of the present step of calculations is to express the last unknown SCALAR R 
through the other SCALARS which may be regarded in this context as known functions. More 
precisely, we shall determine df (R,phi ,phi~) (essentially, the 2-dimensional Laplacian of 
R), cf. line 304. R is then obtained by means of a straightforward integrating (the solving 
of Poisson equation). 

First of all, we deduce the equation determining R from the Einstein-Maxwell equations 
by means of the executing appropriate substitutions of SCALAR VALUES. Specifically, the 
instructions 

input: obtain the SPINOR PART OF EINSTEIN-MAXWELL EQUATIONS; 
input: match substitution rules (12) , (17) , (18) with 
input: the ABOVE EQUATIONS; 

input: renew the ABOVE EQUATIONS: 
input: match substitution rules i 
input: the ABOVE EQUATIONS: 

input: renew the ABOVE EQUATIONS; 
input: match substitution rules 
input: the ABOVE EQUATIONS; 

input: renew the ABOVE EQUATIONS; 
input: match substitution rules 
input: the ABOVE EQUATIONS; 

input: renew the ABOVE EQUATIONS; 

input: match substitution rules (9) , (10) , (13) , (14) with 
input: the ABOVE EQUATIONS; 

input: renew the ABOVE EQUATIONS; 

input: match substitution rules (15) , (16) , (8) , (11) with 
input: the ABOVE EQUATIONS; 

input: renew the ABOVE EQUATIONS; 
input: isolate the ABOVE EQUATIONS; 

entails all the field equations which are not immediately satisfied by virtue of the rela- 
tionships introduced in frames of the current problem specification. The result is two 



15) ,(16) with 



(7), (77) with 



(2) , (3) , (5) , (6) , (22) , (33) , (55) , (66) with 
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EQUATIONS. The first of them, Eq. (1), is in fact trivial. It is satisfied, provided the mean- 
ing of rho is taken into account (cf. the comment following lines 289-291). Concerning 
another one, it may be considered just as the equation determining R. The instructions 

input: resolve Eq. (2) w.r.t. df (R,phi ,phi~) ; 

input: show new EQUATION; 

yield in particular the following the output: 
output: Eq. (3) reads DF(R,phi,phi~) = ( 

output: 4 4 2 

output: - 8*DF(k,phi) *DF(k,phi~) *psi *psi~ *rho 

<g; the continuation dropped out consists of 119 non-empty lines ^> 

(r.h.s. here involves neither R nor its derivatives). The system of field equations in the 
representation described by the current problem specification entails no other relationships 
which involve SCALAR R. 

Although Eq. (3) seems to be fairly involved, it is nevertheless quite tractable. We 
apply the method of the 'partial integrating' which was, essentially, already used above. 
The corresponding instructions are the following: 
input: evaluate 



420 input: expr=R, 

input: aux=LHS (EQ (3) ) -RHS (EQ (3) ) , 

A22input: expr=expr +k**2*C0EFFN (COEFFN (NUMR (aux) 

••• input: ,df (k,phi) ) ,df (k,phi~) ) 

424 input: /DENM (aux) /2 , 

425 input: aux=DR (expr , phi , phi~ ) - (LHS (EQ (3) ) -RHS (EQ (3) ) ) , 
input aux=MATCHING(aux,8) , 

427 input: expr=expr-k*DF (psi~ ,x) 

■■■input: * *C0EFFN (COEFFN (NUMR(aux) ,df (k,phi) ) 

••• input ,df (psi~ ,x, phi") ) /DENM (aux) 

■■■input: -k*DF(psi,x) 

■■■input: *C0EFFN (COEFFN (NUMR(aux) ,df (k,phi~)) 

432 input: ,df (psi ,x,phi) )/DENM(aux) , 

433input: aux=DR(expr , phi , phi~) - (LHS (EQ (3) ) -RHS (EQ (3) ) ) , 

input: aux=MATCHING(aux,8) , 

435 input: expr=expr-k*psi~ 

■■■input: ' *C0EFFN (COEFFN (NUMR(aux) ,df (k,phi) ) 

••• input: ,df (psi~ ,phi~) )/DENM(aux) 

■■■input: -k*psi 

■■■input: *C0EFFN (COEFFN (NUMR(aux) ,df (k,phi~)) 

440 input: , df (psi ,phi) ) /DENM(aux) , 

441 input: aux=DR(expr , phi , phi~) - (LHS (EQ (3) ) -RHS (EQ (3) ) ) , 
input: aux=MATCHING(aux,8) , 

input: aux=MATCHING(aux, 17, 18) , 

input: aux=MATCHING(aux, 15, 16,7,77) , 

445 input: expr=expr-beta*beta~ 

■■■input: 1 *C0EFFN (COEFFN (NUMR(aux) ,df (beta, phi)) 

447 input: , df (beta" , phi") ) /DENM (aux) , 

448 input: aux=DR (expr , phi , phi" ) - (LHS (EQ (3) ) -RHS (EQ (3) ) ) , 
input: aux=MATCHING(aux,8) , 

input: aux=MATCHING(aux, 17, 18) , 

input: aux=MATCHING(aux, 15, 16,7,77) , 

input: aux=MATCHING(aux, 19 , 20) , 
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453 input: expr=expr-beta* (int_xi~+gam) 

■■■input: ' *COEFFN(COEFFN(NUMR(aux) ,df (beta, phi) ) 

••• input: ,df (int_xi~ ,phi~) )/DENM(aux) 

■■■input: -beta~*(int_xi+gam~) 

■■■input: *COEFFN(COEFFN(NUMR(aux) ,df (beta" , phi") ) 

458input: ,df (int_xi ,phi) )/DENM(aux) , 

459input: aux=DR(expr , phi , phi") - (LHS (EQ (3) ) -RHS (EQ (3) ) ) , 

input: aux=MATCHING(aux,8) , 

input: aux=MATCHING(aux, 17, 18) , 

462 input: aux=MATCHING (aux , 15 , 16 , 7 , 77) ; 



Their meaning is quite straightforward. Initially, the working variable expr is endowed 

with the value R (line 420). Further it is modified by means of the adding appropriate 

terms (lines 422-424, 427-432, 435-440, 445-447, 453-458) proportional to k**2, k*DF(psi~ ,x) , 

k*DF(psi,x), etc. The necessary coefficients are determined 'automatically', utilizing the 

value of the another working variable aux. The latter is each time equal modulo Eq. (3) 

to the second order derivative DR( expr, phi, phi") of the current value of expr, (lines 

425, 441, 448, 459). During these transformations the properties (in fact, definitions) of 

various variables expressed in the form of substitution rules (see above the section of 

Substitutions) are taken into account. (These are the instructions which involve the 

macro 'MATCHING'.) In particular, after the execution of the instruction displayed in line 

462, aux still equals the derivative DR ( expr , phi , phi") modulo Eq. (3). Moreover, it is 

easy to see that this relationship just expresses (in a modified but equivalent form) the 

essence of Eq. (3) . 

Ultimately, the instructions 
input: turn on the displaying of negative powers; 
input: evaluate aux, expr factoring out R,k,df (psi ,x) ,df (psi~ ,x) ; 

entails the following output: 



output: The expression 

output: aux 

output: amounts to: 

output: 2 -3 -1 -2 

output: - DF(psi,x) *psi *psi~ *rho 

output: 1 -2 ~2 ~2 

output: *DF(psi ,x) *DF(psi~ ,x) *psi *psi~ *rho 

output: 2 

output: 1 -2 _ 2 

output: + *DF(psi,x)*psi *rho *(gam~ + int_xi) 

output: 4 

output: 2 -1 ~3 _ 2 

output: - DF(psi~,x) *psi *psi~ *rho 

output: 1 -2 -2 

output: + *DF(psi~ ,x)*psi~ *rho *(gam + int_xi~) 

output: 4 

output: 1 -2-1-2 

output: + *DF(psi,x,2)*psi *psi~ *rho 

output: 2 

output: 1 -1-2-21 

output: + *DF(psi~ ,x,2)*psi *psi~ *rho *xi*xi~ ^<~^^0.4 s. 

output: 2 4 
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output: The expression 

output: expr 

output: amounts to: 

output: 12 1 -1 
output: R + *k *psi*psi~ *k*DF(psi,x)*I*psi 

output: 4 4 

output: 1 -11 1 

output: + *k*DF(psi~ ,x)*I*psi~ + k*I*( *psi*beta~ + *psi*gam 

output: 4 4 4 

output: 1 11 

output: + *psi*int_xi~ *beta*psi~ *gam~*psi~ 

output: 4 4 J 4 

output: 1 1 

output: *psi~*int_xi) + *(beta*beta~ + beta*gam 

output: 4 4 

output: + beta*int_xi~ + gam~*beta~ + int_xi*beta~) ?<~*?0.2 s . 



In view of all said above the flat Laplacian of expr equals aux which may be considered as 
a known function. This relation is easily integrable. 

It is evident that the result deduced accomplishes the solving of the problem considered. 



3.5 The solution 

Let us summarize the above results re-casting the formulary input data and computer out- 
put to standard mathematical notations. Namely, we found a generic solution of Einstein- 
Maxwell equations which complies with Eq. (15). Its metric is described by the expansion 
(5) where the tetrad of 1-forms 9 a is defined as follows (cf. [9], Eq. (27.46)): 

e° = p-^d^/V^, e 1 = ¥, e 2 = dx/V2, 

6 3 = (dy-2Im(k, 4> d<f>) + (R + Qy-\\iP\ 2 y 2 )dx)/V2 (28) 

while the non-zero components of the undotted spinor of the electromagnetic field comprise 

= ip and <j> 2 = fnp(-S + ip :<p y). (29) 

Here x,y are the real coordinates, is the complex one (whose complex conjugated coun- 
terpart 4> is considered as independent variable), the other symbols denote their functions. 
Specifically, 

P = 1 + II0I 2 , (30) 
S = £ + {3,^+2iiljk, < p+ikip, 4> , (31) 

R = D-\(\(3\ 2 + m 2 ) 

+2 Re [v + i k (T + 7 + P) ~ P (T + 7))) , (32) 

Q = i Re (^-V, x -^(T + 7 + ^)). (33) 
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Here real D = D(x, (ft, (ft) and k = k(x, <ft, (ft) satisfy the equations 



d,* = -\p- 2 \(r l u?-m 2 

-\p- 2 Re (2 V" 1 ,xx + (tT 1 ) ,x (T + 7)) , 

Kn+\p~ 2 k = -|p->|- 2 Im(^-V,.+^(T + 7 + ^)). (34) 

Complex functions £ = £(x , (ft , (ft) , T = T(x,(ft,(ft) are defined according to the following 
relationships: 

£,*=P" 2 (^"%, T,*=£. (35) 

The functions ift = ift(x,(ft) ^ 0, f3 = f3{x,(ft), 7 = 7(2:, (ft), v = is(x,(ft) are arbitrary. 
They are assumed to be holomorphic with respect to the second (complex) argument (and 
smooth with respect to the first one). 

It is worth mentioning that although the description of the solution presented involves 
some unsolved equations, all they are solved by means of subsequent quadratures. Indeed, 
there are three kinds of equations which, dropping out irrelevant complicating details, look 
as follows: 

(***) (l + zz) 2 d2 ^ Z d f +2f = g(z,z). 

Here g = g(z\ ) Z2) denotes some known function holomorphic with respect to the both 
arguments, i.e. admitting a holomorphic extension to some neighborhood of the 'shell' 
{zi = z 2 } C C 2 . Local solutions / of these equations share the latter property. They can 
be represented in explicit form in terms of quadratures as follows: 

1 1 1 

(*) f = z J dtgitz,?), (**) f = zz J J dt 1 dt 2 g(t 1 z,t 2 z), 
00 







+ Re 



dj(z) 2zj (z) 



dz 1 + zz 



the last formula (local general solution of inhomogeneous Helmholz equation on a sphere) 
being of certain independent interest. Here j = j(z) is an arbitrary holomorphic function 
playing role of 'integration constant'. It is worth noting that the second term in (***) 
represents itself the general local solution of the homogeneous version (g = 0) of Eq. (***) 
and thus (***) is its general local solution in non- homogeneous case. In the formulae 
(*), (**) the 'integration constants' similar to j(z) are dropped out since they had been 
explicitly introduced in the above Eqs. (28-35) describing the solution of the field equations. 
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This is just the origin of the arising of the arbitrary functions v, 7, (3. On the contrary, 
the arbitrary function which arises when integrating Eq. (34) (equivalent to (***)) with 
respect to k was not explicitly introduced, being implicitly involved in (the value of) k. It 
is the fifth arbitrary function (the fourth one is ip). 

Thus, in total, the electrovac solution obtained involves five arbitrary functions of two 
variables holomorphic with respect to one of them. 

Concerning the aspect of the physical interpretation of the field configuration described 
by equations (28-35), one should take into account that a huge 'amount' of degrees of free- 
dom the family of solutions possesses makes fairly difficult, if not impossible, to reveal and 
formulate it in a full generality. Nevertheless one may gain some insight into the problem 
analyzing suitable particular cases (subfamilies) of the field configurations described by 
the above equations. That way, a useful simple example is provided by the following their 
representative (cf. [9], Eq. (27.54)): 

g = dx (D (dy + (Rez/(x,0) - \i/>((/>)\ 2 y 2 ) dx) + d ° : " 



h/#)l 2 (i + l0l 2 ) 2 ' 



d(f) A d(f) 

(Here, with regard to Eqs. (28-35), a minor modification of the gauge and some obvi- 
ous elementary re-definitions were carried out). The arbitrary functions il)((j>),v(x,<f>) are 
holomorphic with respect to 0. 

As opposed to the general case, the meaning of the latter field configuration is im- 
mediately manifest: it represents the generalized plane wave spread against the Bertotti- 
Robinson space-time [13] describing itself non-null homogeneous static electromagnetic 
field. Indeed, the Bertotti-Robinson field configuration arises in the case v — 0,ip — 
constant 7^ 0. On the other hand an appropriate limiting procedure corresponding to the 
nullifying the static electric (magnetic) field yields the pp-w&ve metric. 



4 Summary and discussion 

Summarizing, the paper presented exhibits an example demonstrating the application of 
the specialized computer algebra system GRGec to the problem of the searching for so- 
lutions to the coupled Maxwell and Einstein-Maxwell equations. We investigated their 
general class complying ansatz (15). The latter equation imposes certain restrictions on 
the irreducible curvature spinors which give evidence (basing first of all on important 
particular cases) for the associating the configurations under consideration with radiation 
processes. 

Reassigning all the routine calculations to GRGec, we found a generic solution of the 
problem which turns out to involve five arbitrary functions of two arguments (see subsection 
3.5). In agreement, with a priori expectations its particular case can be interpreted as pp- 
wave spread against a background Bertotti-Robinson space-time, see Refs. [13], [9]. This 
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circumstance partially justifies the working term 'radiative electrovacs' (perhaps too wide 
in the context considered) which was used above for the referring to the class of space-times 
investigated. 

The complete problem specification (GRGec input code) which encodes the formulae 
(28-35) from subsection 3.5 and represents the solution found is displayed in Appendix A. 
Appendix B displays the corresponding GRGec output which proves that we actually deal 
with a solution of the relevant field equations 23 . 

As it has been mentioned, the prevailing aim pursued in the present work is the char- 
acterizing of GRGec system from viewpoint of its suitability for a practical application. 
Additionally, we simultaneously attempted, mostly after indirect fashion, to exhibit possi- 
ble styles of its usage and to demonstrate some of its capabilities. In particular, this is one 
of the purposes the inclusion of section 3 which discusses the processing of the problem in 
the complete form. 

It is instructive to mention that the treatment of the problem considered in subsections 
2.2,2.3 began immediately with consideration of a workable GRGec program (we called 
it problem specification, see subsection 3.1). In particular we did not need, essentially, 
to preliminarily learn GRGec input language and data stuctures to a notable extent, re- 
stricting ourselves to several remarks. At the same time it should be mentioned that, in 
principle, GRGec is a fairly complex system whose complete description occupies hundreds 
of pages. One can estimate therefore a prominent clarity of the organization of GRGec 
language which enables one to comprehend the essence of GRGec input without becoming 
absorbed into specific programming- related issues 24 . The content of Appendix A, where an 
integrated problem specification is exhibited, provides a nice demonstration of the latter 
circumstance. 

To be more specific, one may distinguish, in principle, the two kinds of GRGec code. 
The first one mostly follows a traditional imperative style of wellknown programming 
languages. Typical examples of such a sort scripts are displayed in lines 215-222, 361-366, 
etc. Generally speaking, the instructions contained in these fragments of GRGec script 
realize, in a sense, 'low level' mathematical calculations handling symbolic formulae. Their 
evident characteristic property is, as usually, a high degree of a detailed petty control 
by a user depending on intermediate results required. It should be noted however that 
GRGec provides a rather limited set of facilities intended for such a 'lower level' algebraic 
programming. For example, no loops, conditional and branching operators, subroutines, 
etc., are supported. The point is that, though it might seem surprising, they are superfluous 
here in fact. 

A foundation (and, hopefully, a source of potential advantages) of GRGec lives in 

23 For the sake to provide a performance scale for the timing labels shown in the output proto- 
col (as well as for incidental ones occurred above in the fragments of output listings reproduced), we 
note that the IBM/PC compatible computer with AMD386DX 40 MHz CPU and 8 Mb of RAM was 
used. A fiducial point characterizing its (fairly moderate) performance rate is provided by the time of 
the expanding of (a + b + c) 100 by Reduce (without output). The corresponding script reads 'showi 
time$(A+B+C)**100$showtime;'. It reports about 87 seconds of CPU labour. 

24 The developing of new programs requires more detailed knowledge of course. 
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a distinct area. Essentially, in GRGec the control over majority of calculations is based 

mostly on actions of a maximally high level, approaching the one which might be attributed 

to the relationships and notions characteristic of the application field (the geometry and 

the field theory) itself. As examples of the corresponding programming (a code of the 

'second' kind) one could get the instructions displayed in lines 49-51 or 91,92,93. The series 

of instructions displayed below (having no relation to the problem considered in the present 

work) also serves an illustration of such a 'super-high level' coding: 

input: find UNDOTTED WEYL SPINOR and compare it with sample; 

input: classify the ABOVE STUFF; 

input: calculate factorizing denominators UNDOTTED WEYL INVARIANTS 
input: and write them to a disk file; 

Illuminating the grounds of the approach utilized, one of the cornerstones of a 'knowl- 
edge' implemented in GRGec an d concerning the mathematical and physical theoretical 
issues is the collection of so called data objects which model the basic notions originated 
from the geometry and the field theory. In the present work, we dealt with a fairly moderate 
part of them which included TETRAD, UNDOTTED CONNECTION, UNDOTTED EM SPINOR, and 
some others. A full family of data objects implemented so far is rather numerous (although 
still not exhaustive) and also includes such entities as, for example, the DATA SAVED (to 
a disk storage), ALL the TRACES OF ENERGY-MOMENTA (known in the current point of cal- 
culation), COVARIANT DIFFERENTIALS OF WEYL 2-FORMS, and even DUPLICATED CHARGED 
DIRAC DERIVATIVE OF UNDOTTED DIRAC PHI-SPINOR. These are export names of data ob- 
jects which are used for the referring to them in instructions. (There is also possibilities of 
the access to separate components, if any, of a data object. The corresponding examples 
can be found in lines 83, 84, 86, 87 and others.) In order to determine the value of some data 
object one should apply one of the actions 'ELICIT' (from a DATA included in the problem 
specification), 'CALCULATE' (from the other data objects), 'FIND' (i.e. partially 'ELICIT' 
and partially 'CALCULATE'), 'OBTAIN' (used for equations). Depending on the current state 
of environment, the realization of such actions can be rather sophisticated invoking, in 
particular, (the models of) various relationships originated from the physical theory and 
geometry. Further, given the value of a data object, one may use for its transformation 
various substitution rules applying them either globally (by means of the actions EXCITE 
or LET) or locally (action MATCH). The substitutions with the help of Reduce' package 
COMPACT are realised by the action MATCH COMPACTIFYING. The Reduce' polynomial 
factorizer is invoked by means of the action 'RENEW FACTORIZING'. Alternatively, one may 
'RENEW FACTORIZING NUMERATORS' or 'RENEW FACTORIZING DENOMINATORS'. The Reduce' 
routine performing decomposition of rational functions into simple fractions is invoked by 
means of the instruction 'FRACTIONATE (data_object) WITH RESPECT TO 25 (kernel) 26 ' (or 
'FRACTIONATE CONDENSING. . . ' which additionally re-arranges simple denominators into 
factors-SCALARS), etc. 

Generally speaking, we have listed major part of actions intended for the 'active' han- 
dling of data objects. One sees that, in principle, the controlling facilities supported by 

25 These three 'keywords' may be shortened to 'w.r.t. '. 

26 This term closely corresponds to the notion of a kernel introduced and supported in Reduce. 
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GRGec are in no way intricate while they reveal no potential limitations on the complexity 
of the calculation processes underlaid. It may be stated that GRGec successively follows 
a tendency to ensure a maximal clarity and naturalness of the applied code. 

It is also worth noting that GRGec usually does not advertise about the ways of the 
realizing of the instructions performed. The listing displayed in Appendix B is an instruc- 
tive illustration of the latter circumstance 27 . Indeed, one finds no formulae in the output 
issued by the instructions executed. Nevertheless the main result — the confirmation of the 
satisfaction of the relevant equations — is clearly exhibited. (If one wants to be informed 
in more details, an additional instruction to 'type ALL KNOWN' would give rise to a lot of 
formulae.) 

Similarly, as it was mentioned, the basic mathematical formalism which is used for 
majority of calculations discussed in the present work is the calculus of exterior forms 28 . 
However, this circumstance nowhere explicitely manifests itself (except perhaps of instruc- 
tion 38 and lines 43, 46 of output). This may be interpreted as a particular manifestation 
of the intention to release a user from superfluous details. 

At the same time the implementation of exterior calculus (together with the stan- 
dard tensor methods) constitutes itself a foundation of the library of next level routines 
realizing geometrical and physical theoretical relationships. These two interconnected un- 
derlaid 'strata' of GRGec system are relatively independent from their 'manager' — the 
interpreter. On the other hand the narrow application field of programming system is 
determined mostly by the library of applied routines (and, to a less extent, by the ba- 
sic mathematical tools implemented). Thus, having substituted another applied library 
(which should support an appropriate collection of data objects), one obtains the system 
with another application field which is controlled, essentially, in the same way. Thus the 
approach to the integrating of the constituents of GRGec outlined possesses a flexibility 
which makes GRGec a valuable subject of a practical interest worthing a further working 
up. 

Resuming, it may be stated that GRGec system can be successfully applied for practical 
calculations in the field of gravitation theory, demonstrating a high efficiency and excellent 
convenience in an application. It can be also estimated as a promising base for the further 
development of the efficient tools for doing computer analysis of a wide scope of problems 
in the field of theoretical physics. 
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Appendix A 

Here an example of the problem specification (input code processed by GRGec) is 
displayed. It checks the fulfillment of the relevant field equations by the solution described 
in section 3.5. 



Problem Electrovac Metric. 
Data: 

declare COORDINATES x,y ,phi ,phi~ ; 
declare SCALARS psi , beta, gam, nu(x, phi) , 
rho(phi,phi~) , 

k,xi , int_xi, int_del(x,phi ,phi~) ; 

declare REAL x,y,k,rho,del,int_del; 
declare COMPLEX CONJUGATED phi & phi~, 

psi & psi~,beta & beta", gam & gam~, nu & nu~, 

xi & xi~,int_xi & int_xi~; 

SCALAR VALUE follows: rho=l+phi*phi~/4; 

ABBREVIATIONS comprise 

emS= xi +df (beta, phi) +2i*psi*df (k,phi) +i*k*df (psi ,phi) , 
dd_k=-k/rho**2/2 
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+1/(4 i*rho**2*psi*psi~) 
* (aux-C . C . (aux) , WHERE 

aux=df (psi ,x) /psi+psi* (beta~+int_xi~+gam) ) , 
Q=(l/4)*(aux+C.C. (aux) , WHERE 

aux=-df (psi, x) /psi+psi* (beta~+int_xi~+gam) ) , 

R=int_del 

-(beta*beta~ +k**2*psi*psi~)/4 
+ ( aux+C . C . ( aux ) , WHERE 

aux=nu +i*k*(DF(psi,x)/psi -psi* (beta~+int_xi~+gam) ) 
-beta* (int_xi~+gam) ) /4, 
del=-DR(l/psi,x)*DR(l/psi~,x)/rho**2/2 -xi*xi~/4 
- ( aux+C . C . ( aux ) , WHERE 
aux=2 DR(l/psi,x,2)/psi~ 

+DR(l/psi,x)*(int_xi+gam~))/rho**2/4; 

TETRAD comprises 

component | 0=d phi/psi/rho/sqrt (2) , 
component 1 1=C.C. (component I 0) , 
component I 2=d x/ sqrt (2) , 
component |3= (d y -2 IM(df(k,phi) d phi) 

+(R +Q*y -y**2*psi*psi~/4) d x) /sqrt (2); 

UNDOTTED EM SPINOR comprises 
component I l=psi , 

component | 2=rho*psi* (-emS +y*df (psi ,phi) ) ; 

DOTTED EM SPINOR is HERMITEAN CONJUGATED to UNDOTTED EM SPINOR; 
end of data. 
Substitutions : 

(4) rho -> VAL(rho); 

(7) df(xi,phi~) -> DR(l/psi~,x)/rho**2; 
(77) df(xi~,phi) -> DR(l/psi, x)/rho**2; 

(8) df (k,phi,phi~) -> dd_k; 

(9) df (k,phi,2,phi~) -> DR(dd_k,phi) ; 

(10) df (k,phi,phi~,2) -> DR(dd_k,phi~) ; 

(15) df(int_xi , phi) -> xi; 

(16) df (int_xi~,phi~) -> xi~; 

(17) df(int_xi, phi, phi") -> df(xi,phi~); 

(18) df (int_xi~ ,phi ,phi~) -> df(xi~,phi); 

(19) df(int_xi, phi, 2) -> df(xi,phi); 
(100) df (int_del,phi,phi~) -> del; 

end of substitutions. 

Instructions : 

obtain UNDOTTED MAXWELL EQUATIONS; 

match substitution rule (7) with the ABOVE EQUATIONS; 
obtain EINSTEIN-MAXWELL EQUATIONS; 
match substitution rule (4) 

with the SCALAR PART of EINSTEIN-MAXWELL EQUATIONS; 
match substitution rules (9) , (10) , (17) , (18) 

with the SPINOR PART of EINSTEIN-MAXWELL EQUATIONS; 
renew the ABOVE EQUATIONS; 

match substitution rules (7) , (77) , (8) , (15) , (16) with 

the ABOVE EQUATIONS; 
renew the ABOVE EQUATIONS; 
match substitution rules (4) , (100) with 

the ABOVE EQUATIONS; 
renew and type ALL the EQUATIONS; 
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quit ; 

end of instructions. 

>>slang« 
Run ! 



Appendix B 

The output of the processing by GRGe C the code displayed in Appendix A is given. 
(Some initialization messages including the copy of the input script are dropped out.). 



Total time spent amounts to 4.5 seconds. 

Processing of the problem 'Electrovac Metric' 

COORDINATES are listed below: 

x, y, phi, phi" 
SCALARS dependences are shown below: 
nu: (x,phi) 



gam: 
beta: 



(x,phi) 
(x.phi) 
psi: (x,phi) 
rho: (phi, phi ~) 
int_del: (x, phi, phi") 
int_xi: (x, phi, phi") 
xi : (x, phi, phi") 
k: (x, phi, phi") 
nu" : (x,phi~) 
gam": (x,phi~) 
beta": (x,phi~) 
psi": (x,phi~) 
int_xi~: (x, phi, phi") 



(is added as C.C. 

(is added as C.C. 

(is added as C.C. 

(is added as C.C. 



(is added as C.C. 
(is added as C.C. 



xi": (x, phi, phi") 
No unrealized abbreviations have been specified 
Abbreviations processed are listed below: 

emS, dd_k, Q, R, del 

SCALAR VALUE is shown below: 

phi*phi~ + 4 
rho = . 



of the SCALAR nu) 

of the SCALAR gam) 

of the SCALAR beta) 

of the SCALAR psi) 

of the SCALAR int.xi) 
of the SCALAR xi) 



No unrealized form abbreviations have been specified. . . 
No form abbreviations have been specified. . . 
Total time spent amounts to 6.1 seconds. 

** The instructions given will be executed now ** 
==> obtain UNDOTTED MAXWELL EQUATIONS 

...UNDOTTED MAXWELL EQUATIONS have been obtained ?<~*?0.6 s. 

Total time spent amounts to 6.7 seconds. 

==> match substitution rule (7) with the ABOVE EQUATIONS 

UNDOTTED MAXWELL EQUATIONS have been processed 1*~>10.1 s. 

==> obtain EINSTEIN - MAXWELL EQUATIONS 

...SPINOR PART OF EINSTEIN - MAXWELL EQUATIONS has been obtained J<~»J 6 .4 s. 
...SCALAR PART OF EINSTEIN - MAXWELL EQUATIONS has been obtained J<~*J 6 . 7 s. 
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Total time spent amounts to 13.4 seconds. 

==> match substitution rule (4) with the SCALAR PART of EINSTEIN - MAXWELL 
EQUATIONS 

SCALAR PART OF EINSTEIN - MAXWELL EQUATIONS has been processed M 0.2 s. 

==> match substitution rules (9), (10), (17), (18) with the SPINOR PART of 
EINSTEIN - MAXWELL EQUATIONS 

SPINOR PART OF EINSTEIN - MAXWELL EQUATIONS has been processed l*~*l 2.6 s. 

==> renew the ABOVE EQUATIONS 

SPINOR PART OF EINSTEIN - MAXWELL EQUATIONS has been renewed M 3.3 s. 

==> match substitution rules (7), (77), (8), (15), (16) with the ABOVE 
EQUATIONS 

SPINOR PART OF EINSTEIN - MAXWELL EQUATIONS has been processed M 5.7 s. 
==> renew the ABOVE EQUATIONS 

SPINOR PART OF EINSTEIN - MAXWELL EQUATIONS has been renewed l*~*l 5.9 s. 

==> match substitution rules (4), (100) with the ABOVE EQUATIONS 

SPINOR PART OF EINSTEIN - MAXWELL EQUATIONS has been processed M 6.4 s. 

==> renew and type ALL the EQUATIONS 

— > RENEW ALL the EQUATIONS 

SPINOR PART OF EINSTEIN - MAXWELL EQUATIONS has been renewed M 6.5 s. 
SCALAR PART OF EINSTEIN - MAXWELL EQUATIONS has been renewed M 6.7 s. 
UNDOTTED MAXWELL EQUATIONS have been renewed M 6.7 s. 

— > TYPE ALL the EQUATIONS 

SPINOR PART OF EINSTEIN - MAXWELL EQUATIONS 

is satisfied 
SCALAR PART OF EINSTEIN - MAXWELL EQUATIONS 

is satisfied 
UNDOTTED MAXWELL EQUATIONS 

are satisfied 
==> quit 

Total time spent amounts to 20.3 seconds, garbage collection consumed 
137 ms. 

The DISK file for the copying is being closed done 

Quitting 
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